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ABSTRACT 


According to the Office of the Under Secretary of Defense for Acquisition and 
Technology, the Department of Defense currently does not have sufficient means of 
tracking and/or understanding the magnitude of problems associated with musculoskeletal 
(MS) injuries amongst its civilian and active duty personnel. The Naval Health Research 
Center (NHRC), San Diego, CA, has been conducting its own studies and research with 
regard to MS injuries incurred by Naval Special Warfare (NSW) personnel since October, 
1993. One elite group within this community which was targeted for the research is the 
Navy SEALs. The purpose of this thesis is to determine the impact on operational 
readiness using the measure of effectiveness (MOE): mission dependent man-days lost 
over time as a result of MS injuries incurred by Navy SEALs. Because NHRC has been 
seeking a database sufficient to answer such questions related to operational readiness, 
suggestions for software system design are provided. Results for man-days lost over time 


are presented for two missions, two mission segments, and three separate activities. 
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EXECUTIVE SUMMARY 


According to the Office of the Under Secretary of Defense, the Department of 
Defense currently does not have sufficient means of tracking and/or understanding the 
magnitude of problems associated with musculoskeletal (MS) injuries amongst its civilian 
and active duty personnel. Active duty data sources are even more limited than those of 
civilian data sources. Although each branch of the armed forces regularly collects back 
injury data, there have been indications of under-reporting and problems with 
discrepancies. 

The Naval Health Research Center (NHRC), San Diego, CA, has been conducting 
its own studies and research with regard to MS injuries incurred by Naval Special Warfare 
(NSW) personnel since October, 1993. Because Navy SEALs are an elite group within 
the NSW community and undergo what ts considered the most extensive military and 
physically demanding training in the world, they were targeted for this research by NHRC. 
NHRC has been conducting studies and research of the occurrence of MS injuries incurred 
by NSW personnel, and it has been actively seeking a software system design that would 
generate data sufficient to answer questions related to operational readiness. 

Dr. Keith Prusaczyk, Head of Human Performance Department, NHRC, has been 
aggressively involved in this undertaking. Data collected by a software system temporarily 
in use at Naval Special Warfare Group Two (NSWG-TWO) Medical, Little Creek, VA, 
was released to and provided by Dr. Prusaczyk for use in this thesis. 

The purpose of this thesis 1s to determine the impact on operational readiness 
using the measure of effectiveness (MOE): man-days lost over time (for specified 
missions and mission segments) as a result of MS injuries incurred by Navy SEALs. 
Furthermore, because the data were lacking some of the information necessary to more 
accurately compute the MOE, suggestions for data items to be collected as well as for 
software system design are provided. 

Because the software system which generated the data used in this thesis was 
designed using relational databases, the focal point of the experience tour at NHRC was to 
determine the “link” between output files generated and to extract the data from each of 
those files considered necessary in computing the MOE. This resulted in the creation of a 
new database utilized for computations. 

The total number of Navy SEALs is approximately 2000, of which 900 are 
attached to NSWG-TWO. The database used for this thesis covered a 13 month period 


for NSWG-TWO SEALs. During that period, 343 different SEALs received treatment for 
MS injuries. Before computing the MOE, only the information regarding Navy SEALs 
with MS injuries was extracted. Because some SEALs were given limited duty status 
(1.e., no duty, no upper body, no lower body, no jumping, etc.) and/or received treatment 


for their injuries, two separate files were used in the computations: 


e One file called /imduty composed of all SEALs given limited duty status for 

certain injuries. 

e One file called treatment composed of all SEALs receiving treatment for 

certain injuries. 

Sample missions and mission segments provided in NHRC Technical Report 
Number 95-24 titled “Physical Demands of U. S. Navy Sea-Air-Land (SEAL) 
Operations”, were selected to compute the MOE. Each SEAL mission is composed of 
several mission segments, with each segment consisting of several activities. The MOE 
was computed after determining which particular diagnoses could potentially affect 
specific activities performed in each of the missions and mission segments. 

In order to see the variation from mission(s) to mission segment(s), the MOE was 
adjusted to be activity dependent man-days lost over time as a result of MS injuries 
incurred by Navy SEALs. Three sample activities regularly performed by SEALs were 
selected. Results for man-days lost over time are presented for two missions, two mission 
segments, and the three separate activities. 

The methodology implemented in this thesis is easily adapted to any particular 
population at risk and the suggested software system design can be easily adjusted for use 


by any branch of the armed forces. 
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I. INTRODUCTION 


A. MUSCULOSKELETAL INJURIES IN THE DEPARTMENT OF DEFENSE 


According to the Office of the Under Secretary of Defense for Acquistion and 
Technology, [Ref. 1], musculoskeletal (MS) injuries are a serious topic of concern for the 
Department of Defense (DoD). 


Work-related musculoskeletal disorders (WRMDs) include upper 
extremity, neck, back, and leg conditions. Conditions directly related to 
workspace environment, jobs, and tools account for a significant portion of the 
injuries and illnesses reported to the Department of Labor (DoL). Sprains and 
strains account for approximately 44 percent of the injuries and illnesses resulting 
in lost work time U. S. wide (Bureau of Labor Statisitics, 1992). Forty percent of 
these sprains and strains involved the back and other portions of the trunk. 

Data sources in the DoD are limited. Most databases do no collect the 
data needed to accurately describe the magnitude of WRMSDs. Although back 
injury data is routinely collected by all branches of the armed forces, there are 
indications of under-reporting and disturbing discrepancies.... Data provided by 
the Army Safety Center (ASC) and the DoL Office of Workers Compensation 
(OWSP) are among those data sources reflecting such discrepancies and under- 
reporting. 

Active duty data are even more limited. One of the most reliable data 
sources available, however, is information obtained as a result of retirement of 
personnel for medical reasons. The top five medical retirement diagnoses involve 
musculoskeletal conditions and account for the majority of the retirements. [Ref. 
1] 

The DoD Ergonomics Working Group was in the process of developing a template 
for cross-sectional epidemiologic study to be used by each branch of the armed forces in 
order to describe the magnitude of WRMSD problems in the DoD. According to the 
information obtained from this web site, the process started in 1995. One of the 
advantages of the design of this study was that comparable data would be collected across 
all branches of the armed forces. Furthermore, if the study would determine that “more 
accurate data are a concern for the DoD, funding should be allocated to perform the 
study on a larger scale”. [Ref. 1] 

To date, however, there has been no software system developed to generate the 
type of information necessary to reflect the magnitude of the problem of MS injuries 


incurred by personnel under the cognizance of the DoD. 


B. NAVAL HEALTH RESEARCH CENTER 


Since October, 1993, the Naval Health Research Center, Naval Submarine Base, 
San Diego, CA, has been conducting its own research and studies of exercise- and sports- 
related injuries among particular populations at risk considered to be a concern for the 
United States Navy. Of particular interest was the occurrence of musculoskeletal injuries 
incurred by an elite population within the Naval Special Warfare (NSW) community - the 
Navy SEALs (Sea, Air, and Land). Because Navy SEALs are considered to undergo the 
most extensive and physically demanding military training in the world, MS injuries are 
predominant among them. Because SEALs are an elite group within the Navy and DoD 
and require such a selective screening process (both in terms of physical and mental 
qualifications and performance capabilities), this particular population at risk warrants 
even more special attention. 

Dr. Keith Prusaczyk, Head of Human Performance Division, NHRC, has been 
actively involved in the process of establishing a system of data collection which could be 
utilized by all branches of the armed forces in tracking musculoskeletal injuries incurred by 
Naval personnel. The work completed in this thesis is a result of continued research in 
this area which began during an experience tour at NHRC under the direction of Dr. 
Prusaczyk. 


C. PURPOSE 


Each SEAL mission is composed of several mission segments, with each segment 


consisting of several activities. The purpose of this thesis is two-fold: 


e To determine the impact on operational readiness of man-days lost over time 
(for specified missions, mission segments or activities) as a result of 
musculoskeletal injuries incurred by Navy SEALs. 

e To establish which types of data items should be collected in order to more 
accurately determine man-days lost. 


The data used to accomplish this task were collected by a software system, Medical 
Analytic Data System (MADSYS), which was temporarily in use by Naval Special 
Warfare Group Two (NSWG-TWO) Medical, Naval Amphibious Base, Little Creek, VA. 
To accomplish the goals described above, this thesis provides researchers at NHRC with 
insight and suggestions for the design of a software system more conducive to answering 


such questions of concern. It will also provide a better mechanism for the tracking and 


understanding of the magnitude of problems associated with MS injuries among Naval 
personnel. 

The methodology implemented for determining the impact on operational readiness 
by this particular population at risk may be easily adapted to any such population at nsk. 
Suggestions for the software system design and data items to be collected are also easily 


adjusted for use by any branch of the armed forces within the Department of Defense. 





il. BACKGROUND 


In October, 1993, researchers employed within Special Operations Division 
(SOD), Naval Health Research Center (NHRC), San Diego, CA, submitted a proposal 
titled “Exercise-related Injuries of Special Operations Personnel” to the Naval Medical 
Research and Development Center (NMRDC), Bethesda, MD. 


The purpose of this particular work unit was to decrease the incidence and 
morbidity of musculoskeletal (MS) injuries sustained by Naval Special Warfare (NSW) 
trainees and operators. Study A of the proposal outlined the further development and 
maintenance of a computerized, epidemiologic database to provide researchers, Special 
Warfare personnel, and medical providers with descriptive epidemiologic data regarding 
MS injuries among Special Warfare trainees and operators. The data were to be used to : 
1) determine frequency and distribution of MS injury; 2) determine morbidity due to MS 
injury (1.e. recovery time, acute and long-term sequelae, attrition, cost, operational 
impact); 3) aid in health-care planning; 4) guide and support future research efforts 
regarding risk factors, etiology, diagnosis, treatment, and prevention of MS injury; 5) 
evaluate the efficacy of specific interventions; and 6) serve as the information source for 
the NSW community. [Ref. 2: p. 4] 


Data being collected currently by some of the medical facilities affiliated with the 
NSW community are generated by a software system known as SMART (Sports Medicine 
And Research Team). However, a temporary source of data collection provided by a 
software system, Medical Analytic Data System (MADSYS), was in use at NSWG-TWO 
Medical, Naval Amphibious Base (NAB), Little Creek, VA, from December, 1995 
through February, 1997. These data were released to Dr. Keith Prusaczyk, NHRC, and 


provided for use in this thesis. 


A. DATA PROVIDED 


1. Purpose 


Dr. Prusacyzk was particularly interested in using the data provided by MADSYS 
for two reasons: 1) to determine the adequacy of the data provided in answering some 
extremely important questions such as how MS injuries among Navy SEALs impacts 
operational readiness in terms of down-time; and 2) to show the importance of such data 


collection by attempting to answer that question with the given database. 


\s an exercise physiologist, Dr. Prusacyzk regularly participates in and conducts 
studies regarding sports- and exercise-related injuries, primarily amongst SEALs and men 
enrolled at Basic Underwater Demolition/SEALs (BUD/S). He has continued to provide 
valuable information to the people in charge of the training program at BUD/S. Although 
the data provided by MADSYS were collected at NSWG-TWO, the primary command of 
East coast SEAL teams and units in an operational setting, the procedures implemented 
and information provided by the analysis of these data will lend themselves to the training 
arena as well. The occurrence of MS injures is predominantly high at BUD/S. As a 
result, there exists a great need to reduce the frequency of occurrence not only at 
operational commands, but also at the training level, to reduce the amount of time spent 
retraining individuals who cannot complete BUD/S with the class in which they are 
enrolled. 

The United States Navy SEAL teams today, attached to one of two Naval Special 
Wartare Groups, are derived from the Underwater Demolition Teams (UDTs) that were 
established as a result of World War II. The ongin of the SEALs dates back to 
volunteers organized into Navy Combat Demolition Units (NCDUs) in 1943. Over the 
years the mission of the SEALs, although referred to under different names at various 
points in time, has virtually remained the same: to conduct unconventional warfare, 
counter-guerilla warfare and clandestine operations in both blue and brown water 
environments - the requirement for hydrographic reconnaissance and underwater 
demolition. [Ref. 3] 

SEAL teams undergo what is considered the most selective screening process and 
ultimately the most extensive military training in the world. Selection to this prestigious, 
challenging community is a major accomplishment for enlisted personnel and officers alike. 
Once selected, the largest and most difficult obstacle to overcome is graduation from the 
BUD/S training program at Naval Special Warfare Center (NSWC), Coronado, California. 
[Ref. 3] 

During BUD/S students encounter obstacles that develop and test their stamina, 
leadership and ability to work as a team. The training program is designed to push the 
individual beyond what is considered the normal limits of physical and mental abilities for 
the purpose of keeping only those mentally and physically strong enough to perform under 
the most extenuating circumstances - those which would be present during combat 


situations. 


BUD/S is composed of three training phases: (1) Basic Conditioning (nine weeks), 
(II) Diving (seven weeks), and (III) Land Warfare (nine weeks). Phase I consists of 
physical conditioning in the areas of running, swimming, calisthenics, obstacle courses, 
and small boat seamanship. The first five weeks prepare trainees for week six, better 
known as “hell-week”. During hell-week, students undergo five and one-half days of 
continuous extensive physical training with a maximum of four hours of sleep. During this 
week, the men are pushed beyond limits, not only from extreme mental and physical pain 
and exhaustion from the activities themselves, but also from sleep deprivation. The last 
three weeks of Phase I teaches methods of conducting hydrographic surveys and charts. 
Physical training continues in Phase II but the requirements must be met in half the amount 
of time as in Phase I. Phase II focuses on combat Self Contained Underwater Breathing 
Apparatus (SCUBA), including both open circuit (compressed air) and closed circuit 
(100% oxygen). The emphasis is toward basic combat diving techniques - using 
swimming and diving to transport SEALs from launch point to combat object - a skill that 
separates the Navy SEALs from all other Special Operations Forces. Phase III continues 
with physical training, with time requirements even further reduced, and concentrates on 
land navigation, small unit tactics, patrolling, rapelling, infantry tactics, and military 
explosives. The last five weeks are spent on San Clemente Island where techniques 
acquired are applied to real life situations. [Ref. 3] 

After graduation from BUD/S, graduates undergo basic parachute training at the 
Army Airborne School, Fort Benning, Georgia. Navy corpsmen who complete this 
training also attend two weeks of Special Operations Technical Training at Naval Special 
Warfare Center, Coronado, CA, to receive instruction in diving medicine and medical 
skills and 44 weeks additional training at Fort Bragg, NC. [Ref. 3] 

Upon completion of the above training, both enlisted and officer personnel report 
to one of the Naval Special Warfare Commands as part of a SEAL team where they are 
placed in a probationary status for six months. The Commanding Officer, in accordance 
with Military Personnel Manual (MILPERSMAN) articles 1410380 (enlisted) and 
1420160 (officer), recommends official selection to the Navy SEALs. If selected, enlisted 
personnel are awarded a Naval Enlisted Classification (NEC) Code, officers are given the 
appropriate designator, and both are awarded the “trident” insignia, which is only worn by 
SEALs. 


2. Files Generated by MADSYS 


A software system, the Medical Analytic Data System (MADSYS), was a 
temporary source of data collection for researchers at NHRC to closely examine effects of 
MS injuries on NSW personnel. The data were collected at NSWG-TWO Medical and 
released to Dr. Prusaczyk for use in this thesis. Understanding the construction of the 
database and rearranging the data into a more manageable form for analysis in this thesis 
provided the focal point for research during the experience tour at NHRC. 

The MADSYS software system was written using relational data bases. Such a 
design “links” output files generated by a particular index, in this case, “NDXDEMO”, 
which is an index starting at number 101 and representing an identifiable code (much like a 
social security number) of a particular patient. 

Six output files generated by the system as Database (*.dbf) files are the primary 
files used in this thesis. Fields common to most files are: 

e NDXDEMO - patient identifier. 


e NDXVIST - visit index representing visit number which is unique to 
particular patient with a particular visit. 

e NDXDIAG - diagnosis index indicating a diagnosis number unique to a 
particular patient with a particular diagnosis. 

e NDXTRIM - treatment index representing a treatment number unique to a 
particular patient receiving treatment for a particular diagnosis. 

e NDXDUTY - duty index representing a unique number in which a particular 
patient appears on a limited duty status for a particular diagnosis. 

e ENCDATE - date of encounter of injury . 


All the aforementioned index labels started with 101 and incremented by one. The names 
of the files, a sample output of each file, and descriptors of the fields in each file are given 
in Tables 2.1 - 2.6. 


NOXDEMO [ssn _[wawet _amer [oos [sex [Race Janancn [DATEDEWO 
fot |ooo.a0-0000|o0e own [150670 | 1 [__ifoa01- 

ox funeru-tt1 few [sam Jossrv6s [a a] rlorore6 
os [eaa.a2-2022 foot Iw [peoeer | a] a] sfonor.s6 


BRANCH: 
1 = Navy 4 = Army Nationa! Guard 
2 = Marine Corps 5 = Air Force 
3 = Other 6 = Other 
Table 2.1 Maddemo.dbf Sample File with Descriptors 
























a= 





INJTYPE INJMECH|ACTCATE {ACTSPEC 






INDXVIST | ENCDATE |CLINIC VISTYPE ICDTYPE 
103 1 
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CLINIC: VISTYPE: INJMECH: INJTYPE: ICDTYPE: 

1 = Regular 1 = Initial 1 = Overuse 1 = New 1 = Training/MS 

2 = Physical Therapy 2 = Follow-up 2 = Acute 2 = Recurrent 2 = Training/Other 

3 = Sports Medicine 3 = Other 3=N/A 3 = Chronic 3 = Nontraining/MS 

4 = Athletic Trainer 4 = Unknown 4=N/A 4 = Nontraining/Other 


5 = Other 


1=N/A 5 = O Course 9 = Fastrope 13 = Patrolling 

2 = Running 6 = Weightlifting 10 = Rapelling 14 = Boat Ops 

3 = Swimming 7 = Para Ops/Static 11 = Climbing 15 = Diving 
Calisthenics 8 = Para Ops/FF 12 = Skiing 16 = Other 





Table 2.2 Madvist.dbf Sample File with Descriptors 








ENCDATE 
17-01-96 


DIAGSPEC: 










DIAGCODE: DIAGDESC: 
Represents Standard |ICD9Code Diagnosis Descri 


Table 2.3 Maddiag.dbf Sample File with Descriptors 
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RANK: UNITCATE: OPSTAT: 
1-9=E-1 through E-9 1=ST-2 1 = Platoon 
10 = Warrant 2=ST-4 2 = Operator 
11 - 13 = O-1 to O-3 3=ST8 3 = Support 
14 = 0-4 to 0-6 4=SDVTM-2 4 = Other 
15= N/A 5 = GRU-2 
6 = SBU-2 
7 = Other 
Table 2.4 Madhist.dbf Sample File with Descriptors 
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1 = e-stim 6 = Cryotherapy 11 = Other 

2 = Hydrotherapy 7 =KINCOM 

3 = lonophoresis 8 = Hot Packs 

4 = Ultrasound 9 = Exercises 

5 = Phonophoresis 10 = OMT 
Table 2.5 Madtrtm.dbf Sample File with Descriptors 
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yaa tsg 52 tas to f18.01-96 _]25-01-96 | 







5 = No upper body 9 = No swim w/fins 
2 = No Duty 6 = No lower body 10 = No Lifting 
3 = No Jump 7 = No soft sand 11 = Limit weight bearing 





8 = No swim 


Table 2.6 Madduty.dbf Sample File with Descriptors 


B. DATA EDITING/CORRECTION 


1. Creation of a Single Treatment File 


One drawback in using the output files generated by MADSYS is the fact that 
information about a particular patient has to be viewed separately from file to file. In 
order to conduct the appropriate analysis of these data, it was necessary to retrieve the 
pertinent information from each file and combine it properly into one file, known as the 
treatment file. This required determining the index link which, as stated earlier, is the 
patient identifier, “NDXDEMO”. Three files, Madtrtm.dbf, Madvist.dbf, and 
Maddemo.dbf were imported from Microsoft Excel, Version 5.0, as tables into Paradox 
running from Windows 3.1. A query was then created to build one file with the necessary 
information from the three files. The “join tables” command was used joining the fields 
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NDXDEMO, NDXVIST, NDXTRIM, and NDXDIAG. A sample output of the resulting 


treatment file is given in Appendix A. 
2. Editing the Treatment File 


While checking the accuracy of the data, it was determined that some records 
associated with follow-up visits in the Madvist.dbf file were not correct. When opening 
an existing file for a patient in MADSYS, there is an option to select “2” for “follow-up 
for existing condition”. When selected, a screen appears in which the person entering data 
may select the diagnosis code for which the follow-up visit exists. Once the diagnosis 
code is chosen, it appears, by viewing the screen, that the former existing information for 
that diagnosis (injury type, injury mechanism, etc.) is selected. However, unless the 
person entering the data physically reselects all that information by a click of the mouse in 
those fields, the default selections for the last five fields appear in the Madvist.dbf file as 


follows: 


NOXDEMO [Noxvist leNcDATE [cuNic _[vistyeE _[icotyee _[inutyee _INJMECH [acTCcaTE [ACTSPEC 


30-01-96 | al 


Thus, it appears as though the last five fields were all entered as “1”. 

Approximately 60 percent of all follow-up visits did not reveal the correct 
information. Therefore, the match command in S-Plus was used to match all the follow- 
up visits in the treatment file to their corresponding initial visits in the Madvist.dbf file. 
The incorrect information in the treatment file was then replaced by the correct 
information from the Madvist.dbf file. 


3. Lost Data 


Some of the follow-up visits did not have an initial visit associated with them. It 
was determined previously that the system had failed several times throughout the course 
of the thirteen month period of data collection. Because the system had no automated 
timed-backup saving feature, some of the data were lost. Whenever possible, Brandi 
Schober, who entered all of the data at NSWG-TWO Medical, would re-enter the data 
and/or corroborate what the initial information should be. However, it is unknown 


whether all of the lost data were retrieved. 
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4. Retrieving MS Injury Data Incurred by Navy SEALs 


Because data retrieved from NSWG-TWO Medical were collected for all patients 
and dependents entering clinics at the command, it was necessary to extract only the data 
of interest in this thesis; namely, the MS injuries incurred by those patients who are Navy 
SEALs. A description of ICD-9 Codes [Ref. 4], standard diagnosis codes utilized by all 
major medical facilities and hospitals, together with [Ref. 5] were used to determine which 
types of diagnosis codes affected the MS system. Because the standard ICD-9 Codes 
reflect both MS system diseases/injuries which are extrinsic and intrinsic, it was necessary 
to extract only the extrinsic diagnoses. This 1s because researchers at NHRC are primarily 
concerned with those MS injuries which result from extensive overuse or changes in 
training schedules and exercise, because those are the injuries which may be avoided with 
proper training. Once determined, only those records which reflected extrinsic MS 
injuries were retrieved from the existing treatment file. Of those records, only the patients 
who are Navy SEALs (indicated by an entry of 1 or 2 in the OPSTAT field of the 
Madhist.dbf file) were retrieved. 


5. Changing the Date Format 


The final step in editing the updated treatment file was to change the date format. 
MADSYS generated dates in the format 03-01-96 representing the date January 3, 1996. 
Because neither Microsoft Excel nor S-Plus (the two software packages utilized in this 
thesis) recognized this format as a date format, a macro was written in Microsoft Excel, 
version 5.0, to change the date format to 1/3/96, which is recognized by both software 
packages. A copy of the macro, Function fixdates is given in Table 2.7. 

There was no interest in retrieving the records in the Madduty.dbf file for which a 
patient was given a “full” duty status. (It is important to note that not all patients who 
did not receive a limited duty status received a “full” duty status). Therefore, only the 
records from that file with entries in the DUTYSTAT field not equal to “1” were 
extracted. The new file was named limduty. Both the new treatment and limduty files 
were the two baseline files used to proceed in determining the man-days lost over time as 


a result of MS injuries incurred by Navy SEALs. 


Function fixdate(indate As Variant) As String 
which = Mid(indate, Len(indate) - 2, Len(indate) - 2) 
If which ="-" Then 
fixdate = Mid(indate, 4, 2) & "/" & Mid(indate, 1, 2) & “/' & Mid(indate, 7, 2) 
Else 
myyear = Right(indate, 2) 
inlen = Len(indate) 
If (Mid(indate, inlen - 4, 1) = "/') Then ‘month has one digit 
mymonth = Mid(indate, inlen - 3, 1) 
myday = Left(indate, inlen - 5) 
Else 
mymonth = Mid(indate, inlen - 4, 2) 
my day = Left(indate, inlen - 6) 
End lf 
fixdate = mymonth & "/" & myday & "/’ & myyear 
End lf 
End Function 





Table 2.7 Function fixdates 


C. IMPACT ON OPERATIONAL READINESS 


Musculoskeletal injuries have been identified specifically by the NSW community 
as a high-priority topic for research. During Operation Desert Storm, MS injuries were the 
second most common cause of troop morbidity. At the 1991 and 1992 NSW Biomedical 
Research Conference, exercise-related injuries were listed as the top priority for future 
medical/physiological research. In fact, seven of the eight NSW commands listed 
exercise-related injuries as the most important problem needing attention. In a NHRC 
survey of East (N=88) and West (N=84) Coast SEALs, exercise-related injuries were 
among the top three areas targeted for biomedical research. [Ref. 2: p. 10] 


Obviously, because of the extensive screening and selection process, training and 
subsequent requirements to complete BUD/S and earn the “trident” insignia, the Navy 
SEALs are considered an elite group within the special operations forces and the Navy 
itself. Because so many men may initially be eligible but so few are selected, the Navy is 
obligated to ensure that these highly trained SEALs can perform in the event of combat to 
the fullest of their abilities. 

Fiscal Year 1997 authorization for Navy SEALs was just over 2000 men. It ts 
conceivable, therefore, that man-days lost due to injury of Navy SEALs has a much 
greater impact on operational readiness than loss by other groups because of the small 


population size. 
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D. THESIS ORGANIZATION 


Chapter III describes the process by which steps were taken to compute the crucial 
measure of effectiveness (MOE): man-days lost over time (for specified missions and 
mission segments) as a result of MS injuries incurred by Navy SEALs. The procedure 
entails a new database which was created that represents all the SEALs out on a particular 
day due to limited duty restriction(s), as a result of receiving treatment , or both. A 
description of how MS diagnosis codes were “mapped” to particular limited duty codes 
(such as “no duty”, “no lower body”, “no upper body’) as given in the Madduty.dbf file is 
provided. Discussion continues with selection of specific missions and mission segments, 
and finally, with development of functions to compute the MOE. 

Chapter IV contains the graphical and statistical results obtained by implementing 
the methodology described in Chapter III. 

Chapter V provides NHRC with recommendations and suggestions for future data 
items to be collected in order to more accurately compute the MOE. This chapter also 


provides suggestions for ongoing and future research. 


Hl. METHODOLOGY 


A. DATA COLLECTION ISSUES/ASSUMPTIONS 


Musculoskeletal injury data for Navy SEALs at NSWG-TWO were to be used as 
the data set to determine the impact on the MOE: Man-days lost over time due to MS 
injuries incurred by Navy SEALs. However, it was determined that these data are not 
sufficient to answer the question. To further illustrate the point, recall the data provided 


when a SEAL visits a clinic for a particular injury: 


NDXDEMOINOXVIST [ENCDATE JCLINIC/VISTYPE JICOTYPEJINJTYPE JINJMECHIACTCATE |ACTSPEC |VISTDATE 


Several important pieces of information are missing. For example, there is no 
healing time estimate for the injury. Nowhere is it recorded that the patient is required to 
return for physical therapy for three weeks. The only information provided involving time 
is the actual date of the visit and the date of the last visit for the same diagnosis (not 
shown above). 

There are some instances where patients are given a limited duty status which is 
output to the separate Madduty.dbf file. Recall the pre-defined categories for limited duty 
status that appeared in the Madduty.dbf file generated when the duty status portion was 
completed by the reviewing physician: 


DUTYSTAT: 
1 = Full duty 3 = No upper body 9 = No swim w/fins 


2 = No duty 6 = No lower body 10 = No lifting 
3 = No jump 7 = No soft sand 11 = Limit weight bearing 
4=No dive 8 = No swim 





The limited duty is related to the specific activity that should not be performed. 
However, as seen above, the DUTYSTAT information recorded is not the same as the 


information regarding types of activities in which injury may be sustained: 


ACTSPEC (activity specific): 
5 = O Course 9 = Fastrope 13 = Patrolling 


6 = Weightlifting 10 = Rapelling 14 = Boat Ops 
3 = Swimming 7 = Para Ops/Static 11 = Climbing 15 = Diving 
4 = Calisthenics 12 = Skiing 16 = Other 
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For example, a patient with an ankle sprain may be told by the physician not to run 
or Jump for two weeks, but “no running” does not appear as a duty status category. Thus, 
some consideration was given to make the activity and limited duty status categories 
reflect the same information, which could also easily be tracked to each of the missions 
and mission segments used during the analysis. As a result, a twelfth DUTYSTAT 
category was created: “12” for no running/walking with additional stress factors (such as 
carrying a 100 lb pack while trying to walk). This category was necessary because almost 
all of the missions and mission segments performed by SEALs result in walking under very 
stressful conditions. Although running seldom appears in the description of missions and 
mission segments (provided for use in this thesis) performed by SEALs, it is conceivable 
that this activity may be necessary when performing missions during combat operations. 
For purposes of this thesis, the other types of training activities which should not be 
performed (such as rapelling, fastroping, etc.), were not given limited duty status codes, 
since the database is a reflection of SEALs injured while in an operational status. 
However, this type of information is easily added to reflect any environment impacting the 
MOE. 

A patient given limited duty status for two weeks may have continued to see a 
physical therapist every day for three months. This implies that although no longer limited 
in duty, the healing process actually was not complete. Clearly, the individual expected to 
engage in these activities during combat would not have the luxury of being treated by 
physical therapy on a daily basis and may not be able to perform this activity (or perform 
the activity to his maximum ability). This would certainly affect his ability to perform any 
of the missions or mission segments which require use of the particular activity. For 
purposes of this thesis, it was assumed that an individual receiving treatment for a time 
period, but not given a limited duty status, would only be able to perform the activity at 50 
percent of his normal ability, while those on limited duty would not be able to perform the 
activity at all. The percent effectiveness of individual performance can easily be changed 
by anyone using the database to compute the MOE. 

It is unknown whether the last visit for a patient with a particular diagnosis 
indicates the patient has actually been healed, is no longer part of the command, or if that 
was the last visit for that diagnosis before data collection ceased. Because this 
information was not available, it was assumed that the patient’s duration of treatment for 
particular injury started at the time of his diagnosis and ended with his last visit for that 


diagnosis. 
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B. CREATION OF TREATMENT/LIMITED DUTY DATABASE 


1. Treatment and Limited Duty Files 


Two files which already existed, the sreatment and limduty files, were modified for 
ease and efficiency of computation. For instance, the treatment file showed every visit for 
which a patient received treatment for each diagnosis, resulting in the same information 
being repeated, with the exception of the date of visit. As mentioned earlier, the 
assumption was that the duration of treatment was from the time of diagnosis for injury 
through the last treatment for the diagnosis. A function, diag.days (see Table 3.2), was 
written in S-Plus to query the treatment file for “START” and “STOP” dates of 
treatment, and to call another function, todate (see Table 3.1) to change the date format 


function(mmddyy.vec) { 
# Convert a date like 1/13/95 into a number like 950113. 
result <- numeric(length(mmddyy.vec)) 
for(i in 1:length(mmddyy.vec)) { 
mmddyy <- mmddyy. vec[i} 
ne <- nchar(mmddyy) 
yy <- as.numeric(substring(mmddyy, ne - 1, nc)) 
if(substring(mmddyy, ne - 4, ne - 4) == "/") { 
# one-digit day dd <- as.numeric(substning(mmddyy, ne - 3, nc - 3)) 
if(ne > 6) 
mm <- as.numeric(substnng(mmddyy, 1, 2)) 
else mm <- as.numenc(substnng(mmddyy, 1, 1)) 


dd <- as.numeric(substring(mmddyy, nc - 4, nc - 3)) 
if(ne > 7) 


mm <- as.numenc(substnng(mmddyy, 1, 2)) 
else mm <- as.numenc(substnng(mmddyy, 1, 1)) 
} 
result[i] <- yy * 10000 + mm * 100 + dd 
} 


return(result) 





Table 3.1 Function todate used to change date format to “yymmdd” 


from “mm/dd/yy” to “yymmdd”. All index categories and diagnosis codes were retrieved. 
The resulting file was named tridays (see Appendix B). 
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function(data = treatment) 
{ 

cols <- c¢(/NDXDEMO.2", "NDXDIAG.1", "DIAGCODE”, "First', "FirstSre", "Last’, "LastSrc", 
Diff) 

result <- matrix(0, length(unique(data[, "NDXDIAG.1"])), length(cols)) 

dimnames(result) <- list(NULL, cols) 

result <- as. data frame(result) 

u.diag <- unique(data[, "NDXDIAG. 1"}) 

for(i in 1:length(u.diag)) { 
cat("Working on diagnosis ", u.diag[i], "\n") 
rows <- (1:nrow(data))[data[, "NDXDIAG.1") == u.diag[i]] 
first.row <- data{rows[1], , drop = F] 
vec <- c("NDXDEMO.2", "NDXDIAG.1") 
result{i, vec] <- first.row[vec] 
result[i, "DIAGCODE"] <- as.character(as. vector(first.row['DIAGCODE"]$DIAGCODE)) 
e.date <- todate(as.character(data[rows, "ENCDATE"))) 
d.date <- todate(as.character(data[rows, "DIAGDATE"))) 
overall.min <- min(e.date, d.date) 
overall.max <- max(e.date, d.date) 
result{i, "First"] <- overall.min 
result{i, “Last'] <- overall.max 
if(min(e.date) < min(d.date)) 

result[i, "FirstSrc"] <- "Enc" 
else result{i, "FirstSrc”] <- "Diag" 
if(max(e.date) > max(d.date)) 
result{i, "LastSrc"] <- "Enc" 

else result{i, "LastSrc"] <- "Diag" 
min.yy <- trunc(overall.min/10000) 
max.yy <- trunc(overall.max/10000) 
min.mm <- trunc((overall.min - 10000 * min.yy)/100) 
max.mm <- trunc((overall.max - 10000 * max.yy)/100) 
min.dd <- overall.min - 100 * as.integer(overall.min/100) 
max.dd <- overall.max - 100 * as.integer(overall.max/100) 
yys <- max.yy - min.yy 
mms <- max.mm - min.mm 
dds <- max.dd - min.dd 


if(mms < 0) { 
mms <-mms + 12 
yys <- yys-1 

} 

if(dds < 0) { 


mms <- mms - 1 
dds <- dds + 30 


} 
result{i, "Diff'] <- paste(format(c("12", yys)[2]), "/', format(c("12", mms)[2]), 
/, format(c("12", dds)[2]), sep = "") 


} 


return(result) 


Table 3.2 Function diag.days used to create a treatment file with start and stop dates 


The /imduty file already had fields containing the “START” and “STOP” dates of 
limited duty, but did not record the diagnosis code. This additional information was 
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appended to the /imduty file using the match function in S-Plus, which matched each 
diagnosis index (NDXDIAG) in the dimduty file to the same diagnosis index in the ¢rtdays 
file (since the NDXDIAG 1s unique to each patient). Sample output for each of the files 
can be seen in Tables 3.3 and 3.4. Note that the last column of the ¢ridays file, “Diff’, 


represents the total time in treatment in “years/months/days’’. 





_ 102] 981726691000 | 961102|Diag | 961107|Enc__|o/o/5 


Table 3.3 Sample output of fridays file 







NOXDEMO_ INDXVIST [NDXDIAG INDXDUTY [OUTYSTAT [START {STOP _{DIAGCODE _ 
ist 1126 s*L'951204__ |960109 924200000 __| 


oe _| 
za [ 2821 “fuse o aso 13 foe0121 700000 — 
2g (133129124 Ci 960712 _]960126 726691000 __ 
Table 3.4 Sample output of limduty i 





2. Database Design 


The design of the database to be used for computing the MOE must show, for 
each day for which data were collected, the individuals who were in treatment and/or on 
limited duty for a particular diagnosis. The dates of data collection were from December 
4, 1995 through February 11, 1997 (951204 - 970211). 

There is also the need to separate the treatment information from the limited duty 
information when performing computations. Therefore, three pertinent pieces of 
information are retrieved from each of the fridays and limduty files: 


e Patient identifier, Demo. 

e Category, Cat, “T” representing the data were retrieved from the frtdays file or 
“D” indicating retrieval from the /imduty file. 

e Diagnosis code, Diagcode. 


A function make.table (see Table 3.5) was written in S-Plus to retrieve, for each 
day, the three field records first from both the trtdays and limduty files. The table was 
created by a call to another function, out (see Table 3.6) also wnitten in S-Plus. The 
primary purpose of the make.table function is to keep track of the days of the month in the 


database and rely on “out” to “make the table”. The result was a non-rectangular database 
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consisting of each row representing a day. Each row contained the new date in “yymmdd” 
format, followed by the three field records for all the “T” categories, followed by the 
three field records for all the “D” categories. 


function() 


{ 


# Calls function "out" to look for patients receiving treatment and/or on limited duty status for each day 


# "Out" passes back a row vector for each day containing the date and 3 field records 
# Make table of IDS of SEALs out because of their being 
# on actiuve treatment or limited duty, by date, two 
# columns per day of 1996. 
today <- 951204 
moday <- data.frame(Mo = 1:12, Days = ¢(31, 29, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31)) 
yy <-95 
mm <- 12 
dd <- 4 
i <- 1 
repeat { 
cat("Date is now ", today, “\n") 
cat(out(today)) 
if(today >= 970211) { 
break 
} 
i<-i+1 
today <- today + 1 
dd <- dd + 1 
if(dd > moday[moday[, "Mo"] == mm, "Days’)) { 
dd <- 1 
mm <- mm + 1 
if(mm > 12) { 
yy <-yy +1 
mm <- 1 
} 
today <- 10000 * yy + 100* mm +dd 





Table 3.5 Function make.table used to create the dutyntrt file 
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function(day) 
{ 
# Out: Which SEALs are out on a given date? 
# argument: day, the relevant date in yymmdd format. 
# This function looks at the two tables "limduty” and 
# "trtdays" and produces a vector of |D numbers 
# (from NDXDEMO) of people out on this day. 
trtvec <- tridays[, "First"] <= day & trtdays[, "Last"] >= day 
badt.ttt <- trtdays|trtvec, "NDXDEMO"] 
diagt.ttt <- as.character(as. vvector(trtdays|trtvec, "DIAGCODE"))) 
Idutyvec <- limduty[, "START"] <= day & limduty[, "STOP"] >= day 
badt.duty <- limduty[Idutyvec, "NDXDEMO"] 
diagt.duty <- as.character(as.vector(limduty|[Idutyvec, "DIAGCODE"})) 
if(sum(trtvec) == 0) 
ttvec <- NULL 
else ttvec <- paste(",T", badt.ttt, diagt.ttt, sep = ",", collapse = “") 
if(sum(Idutyvec) == 0) 
Idvec <- NULL 
else Idvec <- paste(",D", badt.duty, diagt.duty, sep = ",", collapse = "") 
paste(day, ttvec, Idvec, "\n", sep = '") 





Table 3.6 Function out used to obtain data on SEALs out each day 


The resulting database in S-Plus was “dumped” from S-Plus, because it was not a 
rectangular table (since the number of patients in treatment and/or on limited duty status 1s 
not necessarily the same from day to day). This was done using the following commands 


at the prompt: 





> sink ("c:/snowth~1/dutyntrt") # allocates space for file in snowthesis directory 
> make.table () # creates file 
> sink # copies the file to snowthesis directo 










These commands allow the user to open the file (dutynirt) created using any 
applicable software package. This particular file was opened in Microsoft Excel. Because 
the total length of each row was unknown, no category headings were assigned prior to 
opening the file in Excel. Category headings were inserted into the file in Excel. In order 
to perform future computations in S-Plus, it was necessary to retrieve the file from Excel. 
However, S-Plus does not recognize non-rectangular tables. Therefore “dummy” records 
were inserted in order to make the file readable in S-Plus. The “dummy” records would 


prove useful later (in terms of efficiency and run-time) in searching the dutynirt file for 


zt 


individual records. This is because during the search, once a dummy record was 
encountered on a particular row, that row would be exited and the search would proceed 
to the next row. The updated dutyntrt file was read into S-Plus using the read.table 
command. Sample output of the created database, dutyntrt, is given in Table 3.7. A 


dummy records appears in the last three fields of row one in the sample output. 





Day _|Cat|Demo |Diagcode _|Cat |Demo |Diagcode _|Cat|Demo |Diagcode _|Cat|Demo |Diagcode_| 
'951231|T_ | 115|924200000|T | 252|8310A0000|D | 115] 924200000/N_| 999] 99999) 
'960101|T_ | 115|924200000|/T_ | 252|8310A0000|D | 115] 924200000/D | 215] 822000000 
960102|T_ | 115}924200000[T_| 212] 844200000|T_ | 25218310A0000 |p | __115] 924200000] 


Table 3.7 Sample output of dutynitrt file 







C. MAPPING DIAGNOSIS CODES TO LIMITED DUTY (LD) CODES 


In order to determine how the SEALs with particular diagnoses on a given day 
affect performance of specific missions and mission segments, it was first necessary to 
determine how each of the diagnosis codes affects an individual’s ability to perform the 
specific activities involved. These activities include those for which the individual is given 
limited duty status (jumping, diving, upper body, lower body, etc.). Of course, the manner 
in which a certain injury affects an individual varies with each individual. In fact, it is 
possible that one person with a particular injury would still be able to perform an activity 
with limited capability while another may not be able to perform the activity at all. 
However, since this type of information - specific to each individual sustaining injury - 
was not available, a more subjective approach was taken in order to determine how each 
of the diagnosis codes could affect each of the LD codes. 

First, all the unique diagnosis codes were extracted from the treatment and limduty 
files. They were then imported into a file, wnqcode.x/s, in Microsoft Excel. Each row 
label was a diagnosis code (in ascending order) and each column label was “LD2” through 
“LD12” representing the limited duty status categories “2” through “12”. Previous 
information obtained from patients given limited duty status was used to assist in the 
determination of which diagnoses affect which activities (e.g, a majority of the back 
injuries resulted in “no lower body”). Therefore, the assumption was made that those 
back injuries which resulted in a patient being unable to use his lower body could (in a 
worst case scenario) affect any individual with the same back injury. Keep in mind that it 
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is possible that any of the injuries could result in “no duty’, “no lower body” or “no upper 
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body”. However, in order to obtain a more realistic mapping from diagnosis to LD code, 
only those injuries which were the most severe, such as injuries associated with the 
vertebral column, were assumed to result in “no duty” (LD code 2). The same approach 
was used to determine which diagnoses could result in “no upper body” (LD code 5) and 
“no lower body” (LD code 6). These three limited duty status categories were of 
particular importance because they are the ones which result in the largest fraction of man- 
days lost over time due to MS injuries. 

In cases where no previous information existed to aid in determining the mapping, 
both [Ref. 5] and Brandi Schober (physical therapist who entered all data at NSWG- 
TWO) provided valuable insight and information as to which types of injuries could 
potentially result in which types of limited duty activities. Each diagnosis was investigated 
in [Ref. 5] to determine the exact muscles and joints affected by the MS injury. 

One particular code, “V67000000”, referring to post surgical procedure does not 
indicate for which injury the surgery was performed. Therefore, it is necessary to extract 
the patient records from the /imduty and treatment files where this code appears and 
compare the information to that in the Maddiag.dbf and Madvist.dbf files to replace that 
code with the code for which the surgery was required. Any records in the files for which 
these codes resulted in a non-musculoskeletal injury were removed from the two files. 

Once the determination was made as to which diagnosis codes mapped to 
particular LD codes, the next step was to reflect this information in the ungcode file. This 
was done by assigning a “1” in the column of the LD codes in which the corresponding 
row entry, diagnosis code, could affect it; otherwise, a zero was assigned. This is, of 
course, comparable to a Bernoulli distribution where success (1) indicates the LD code is 
affected by the diagnosis code, or failure (0) which indicates no effect. It is important to 
note that this information was determined in the best manner possible, given the current 
database. Any determination of an injury potentially resulting in an LD code which ts not 
reflected in the file can easily be changed by the user. 

The resulting unqcode.txt file was read into S-Plus and saved as object unqcode 
(see Appendix C). A sample output of the file is shown in Table 3.8. 







| 729000000|_o| 1] of of 1] 4] of of o| of 14 
715970000] _o| 1{ of of of 4| of 4] of of 14 
719410000! _o] o| 1] of of of 1] 4] 4] of oO 


Table 3.8 Sample output of unqgcode. txt file 
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D. CREATION OF TABLES FOR EACH LD CODE 


The next step was to write a function in S-Plus which showed all the patients who 
would be affected by each of the LD codes on a given day. Again, the assumption was 
that each diagnosis would result in either the inability or limited ability of the individual to 
perform those activities that the particular diagnosis would affect. For example, as seen in 
Table 3.8, the diagnosis code “729000000” could affect LD codes 3, 6, 7, and 12, which 
means it could result in “no jumping”, “no lower body”, “no soft sand”, or “no running” 
(or a limited ability to perform these). It was necessary to determine which patients were 
put on limited duty and/or received treatment for a particular injury for each LD code (i.e, 
which patients were affected by “no jumping” as a result of being put on limited duty 
status or by “limited ability to jump” as a result of receiving treatment). 

A separate function was written for each LD code which would search the dutynirt 
file each day. This was done not only to query the dutynirt file for one LD code at a time, 
allowing for timeliness in search, but also to allow for ease of editing records in any one 
table. For example, suppose it was discovered that an error had been made in column 
three of the ungcode file. Only those patients affected by LD code 3 would change. A 
separate function for each LD code allows the user to correct the information in the 
ungqcode file and rerun the function without having to change any of the other tables 
already created. For every three field record (Demo, Cat, Diagcode), the function 
searches the ungcode file for that particular diagnosis code and outputs a four field record 
(Day, Demo, Cat, Diagcode) for each patient affected by that particular LD code. For 
example, function /dtable2 searches the dutynirt file for all the patients affected by LD 
code 2 and outputs the four field record into a table called table2. Each function was 
named /dtableX, where X represents the appropriate LD code number. Each function 
would return tableX ,where X is as defined above. Each of the tables returned are given 
column headings and saved as /drX.out (LD results), where X is defined above. Function 
Idtable2 is given in Table 3.10. A sample output of the table /dr2. out is given in Table 3.9. 







960220 [D_|231__|721900000_ 


Table 3.9 Sample output of /dr2.out 
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Next, each of the /drX.out files are “filtered” in S-Plus to obtain all those with 
category “I” separate from all those with category “D”. Each is saved as /drXt. out (those 
on limited duty X) or /drXd. out (those affected by limited duty X but receiving treatment). 
Before any computations of man-days lost over time for any of the missions and mission 
segments are made, it 1s necessary to determine which missions and mission segments are 


required. 


function(rowcodes, unqcode) 
{ 
# rowcodes here is the dutyntrt file 
# unqcode is the unqcode file 
t2.ctr <- 1 
table2 <- matrix(" ", nrow = 2500, ncol = 4) 
rowcodes.col <- ncol(rowcodes) 
for(i in 1:nrow(rowcodes)) { 
cat("Now on row “, {) 
j.seq <- seq(4, rowcodes.col, by = 3) 
js <- rowcodes{i, j.seq] 
js <- jsfjs = "99999" & js !=" 99999"] 
j.seq <- seq(4, by = 3, length = length(js)) 
samprows <- unqcode{[as.logical(match(unqcodef, 1], js, 0)), , drop = F] 
which.are.2 <- samprows[samprows{, 2] == "1", 1] 
if(length(which.are.2) != 0) { 
cat("...found!\n") 


j.seq <- j.seq{as.logical(match(js, which.are.2, 0))] 


right <- cbind(i, c(rbind(1, j.seq - 2, j.seq - 1, j.seq))) 
n <- length(j.seq) 
left <- cbind(rep(t2.ctr:(t2.ctr + n - 1), rep(4, n)), rep(1:4, n)) 
table2{left] <- rowcodesj{right] 
t2.ctr <- t2.ctr +n 
} 
else cat("...nothing there\n’) 


} 
return(list(table2 = table2)) 





Table 3.10 Function /dtable2 to show all patients affected by LD code two 
E. SEAL MISSIONS AND MISSION SEGMENTS 


As previously mentioned, the next step in computing the MOE is to select a sample 
of several missions and mission segments performed by Navy SEALs. A study 
documented in [Ref 6] sought to determine, by interviewing 82 SEALs from NSWG-ONE 
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and NSWG-TWO having an average time in community of approximately 11 years, what 


types of activities and abilities were considered by SEALs to be the most crucial in 


performing certain missions and mission segments. The primary missions were 


categorized broadly as: 


Beach surveys in support of amphibious operations. 

Reconnaissance and intelligence gathering. 

Rescue of downed pilots and extraction of personnel from dented territory. 
Limpet (ordnance) attacks against enemy ships and piers. 

Direct actions against military targets on shore or at sea. 


The missions are composed of several segments which include planning, transit, 


insertion, Operation, extraction, and egress. Each segment requires specific, distinct 


activities that are integral components of the overall mission. 


A survey questionnaire was answered by the participating SEALs. The survey 


questions concerned three areas of description: 


Missions. 
Mission segments. 
Thermal stress conditions (not used 1n this thesis). 


Section (1) described 37 missions falling into five categories: 


Overland - 13 of primarily walking and carrying heavy loads. 

Across the Beach - 11 involving transit from a cold, wet environment to a cold, 
dry environment. 

Ship Attack - 6 limpet attacks incorporating varied insertion modes. 

Survey - 2 beach surveys. 

Other - 5 encompassing attack on an oil platform, equipment loading onto a 
submarine, and 3 combat missions for which SEALs received Medals of 
Honor. 


Section (2) described 50 mission segments falling into four categories: 


Walking/hiking/skiing - 12 of primarily walking carrying heavy loads. 
Swimming/diving - 17 surface or underwater swimming using SCUBA. 
Lifting/dragging/carrying/climbing - 16 involving rescue carriers and carrying 
heavy equipment and climbing. 

Jumping/bumping - 5 using fast rope or parachute for insertion. 


SEALs were asked three questions about each section and were instructed to 


record a number from a seven point Likert scale (values 1 - 7) that best reflected their 
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opinion and personal experiences. Example scale levels were of the form “easiest” (=1), 
“medium difficulty” (=4), and “most difficult” (=7) with intermediate values allowed. In 
addition, SEALs were asked about the frequency with which such missions were 
performed: “once per career’, “once per year”, or “once per month (or more)”. SEALs 
were also asked to rank them in order of importance. Mean scores were computed for 
each of the categories (difficulty, frequency, and importance). Each mission and mission 
segment was ordered in terms of difficulty, frequency, and performance. An “overall 
criticality”, the sum of the mean scores for each of the categories, was computed. All 
missions and mission segments as categorized above are given in Appendices B through I 
of [Ref. 6]. 

A few missions and mission segments which ranked high in two of the three 
categories (difficulty, frequency, importance) were selected for use in this paper. Also, the 
selection was made to result in some of the same activities, as well as some distinct 
activities, being used across the sample. The mission ranked highest in overall criticality, 
(referred to a Mission 1), was used in the sample because it ranked first in frequency, fifth 
in importance and was not considered very difficult to perform (ranked 31 of 37). The 
mission segment ranked highest in overall criticality (referred to as Mission Segment 1) 
was also a good candidate, since it did not use the same activities as Mission 1 and ranked 
first and eighth in frequency and difficulty, respectively. Mission Segment 3 (third in 
overall critically) was selected because it overlapped the swimming portion of Mission 
Segment | and ranked second and fifth in frequency and importance, respectively. One 
other mission was selected that required use of different activities from the previous 
mission and mission segments. As a result, Mission 25, although not ranked high in any of 
the three categories, was chosen primarily because it utilized different activities. A 
description of each mission and mission segment chosen is provided in Table 3.11. The 
percentages in parentheses at the end of each description represent the percentage of 


SEALs interviewed who actually performed the mission or mission segment. 


D4 jj 












Mission 1: 





Walk (hump) 15 km (9 miles) over uneven terrain at night carrying 125 pound pack (with 
radios, etc.) in 70° F temperature to objective; then, retrace steps to extraction point. (85 %) 


Mission 25: 





Parachute (static line) into desert with 150 pound pack (day temp 110° F); then walk 
75 km (47 miles) during the next 7 nights. (25 %) 





Mission Segment 1: 
Swim for 3 hours underwater in 55° F temperature, wearing wet suit, mask fins, and using 
Drager UBA and attack board. (90 %) 





Mission Segment 3: 
Swim/dive a distance of 2000 meters in 56° F water wearing wet suit, fins, and a Drager 
and carrying a limpet mine; then swim 2000 meters without the limpet mine. (89 %) 

Table 3.11 Sample missions and mission segments used in analysis 


F. COMPUTING THE MOE FOR EACH MISSION AND MISSION SEGMENT 


Careful attention was required to assure that man-days lost were not double 


counted in the following situations: 


e A patient receiving both a limited duty status and treatment for the same 
diagnosis on a given day. 

e A patient receiving a limited duty status and/or treatment for one or more 
diagnoses on a given day. 

e A patient affected by more than one LD code for any diagnosis on a given 
day (which is the case for every patient in the database, since every diagnosis 
Is mapped to more than one LD code). 


The easiest way to avoid recounting was to make pairwise comparisons, first 
between each of the /dtables affecting each particular mission or mission segment, then 


between those records involving patients out because of limited duty or treatment. 
1. Comparison of LD Codes for Each Mission and Mission Segment 


Mission 1 is used to describe the process steps in detail. It was determined that 
Mission 1 is affected by LD codes 2, 5, 6, 10, 11, and 12. As a result, the files /dr2t.out, 
ldrSt.out, ldr6t.out, ldr10t.out, ldr11t.out, and ldr12t.out are compared to each other to 
assure no duplicate counting in any of the “T” (treatment) categories. The comparable 
files /drXd.out (X for 2, 5, 6, 10, 11, and 12) are compared to assure no duplicate 
counting in any of the “D” (duty) categories. Once all the pairwise comparisons [ (2,5), 
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(2,6), (2,10), (2,11), (2,12), (5,6), (5,10), (5,11), (5,12), (6,10), (6,11), (6,12), (10,11), 
(10,12), and (11,12) ] are made for all the /drXt.out and ldrXd.out files, the most recent 
version of each of those are compared to each other to ensure an entry that appeared in 
both the “treatment” and “duty” tables for the same day would not be counted twice. 
Functions in S-Plus were written to make each of the pairwise comparisons separately and 
will be discussed later. 

The first level of comparison [ (2,5), (2,6), (2,10), (2,11), and (2,12) ] required 
reading each of the two files to be compared and the addition of another column which 
indicated the LD code. For example, in comparing /dr2t.out and /drSt.out, both files are 
read into function m/t/ (for Mission 1 Treatment Step 1) and stored as /dr2t and /dr5t, 
respectively. Columns of “2’s” and “5’s”, indicating that each record would be accounted 
for by LD code “2” and LD code “5”, respectively, are added to the files. Since LD code 
2 (no duty) would limit an individual more than any of the other codes, the /dr2t file 1s 
used as the main file (i.e., meaning it would not change, but other files may change). Each 
day present in the dutynirt file is “looked for” in each of the two files. A comparison 1s 
made between all records in the /dr2t file and the /dr5t file for a given day. If any entry 
appears in both, then the entry “5” in column 5 of the /dr5t file was replaced by a “2” in 
order to signify it was already accounted for by the /dr2t file. Function m/t/ 1s given in 
Table 3.12. Similar techniques are employed for each of the “first level” pairwise 


comparisons. 


function(Idr2t.out, Idr5t. out) 
{ 


# mission1 t1: deluxe version 


# compares the 2 files and places a "2" 


# in column 5 if a record in Idr5t.out exactly matches one 
# in Idr2t.out 


matchers <- match(vec.5, vec.2) 

Idr2t <- cbind(Idr2t.out, rep("2", nrow(Idr2t.out))) 
Idr5t <- cbind(Idr5t.out, rep("S", nrow(IdrSt.out))) 
new.col <- ncol(Idr5t) 

IdrSt[!is.na(matchers), new.col] <- "2" 
return(list(ldr2t = Idr2t, Idr5t = IdrSt)) 





Table 3.12 Function m/¢/ to compare “treatment” LD codes 2 and 5 
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Second level pairwise comparisons [ (5,6), (5,10), (5,11), and (5,12) ] were then 
made. For example, the function m/t6 (Mission | Treatment Step 6) compares the most 
recent versions of the /dr5t.out (mJoutSt) and the /dr6t.out (m2out6t) files. The most 
recent versions of each of those files now have five columns with either “2” or “5” in 
column five of the m/out5t file and either “2” or “6” in column five of the m2out6t file. 
Again, for each day, a comparison of records 1s made. If the same record (fields one 
through four) appears in both files, then the information is updated in the m2out6t file. In 
other words, if column five of the m2out6t file does not already contain a “2”, then the “6” 
that is currently there is replaced by a “5S” to indicate that it is already accounted for by the 


mloutdt file. Function m/t6 can be seen in Table 3.13. 


function(m1out5St, m2outét) 
{ 
# mission 1 t6: deluxe version 
# Receives the first updated version of the file Idr5t.out (m1out5t) 
# and compares it to most recent version of Idr6t.out (m2out6t); can 
# only use apply on cols 1 thru 4 of m1out5St and m2out6t; places a "5S" 
# col 5 of m2outét if that row exactly matches a row in m1tout5t 

vec.5 <- apply(miout5t[, c(1:4)], 1, paste, collapse = “") 


matchers <- match(vec.6, vec.5) 

Idr6t <- m2out6t 

new.col <- ncol(idrét) 

not2s <- Idrét[!is.na(matchers),new.col] != 2 
# take the entries that haven't already been accounted 
# for in Id2 file and assign them a "5" in col 5 
Idr6t[(!is.na(matchers)) & (not2s), new.col] <- "5S" 
return(list(Idr6t = Idr6t)) 





Table 3.13 Function m/t6 to compare “treatment” LD codes 5 and 6 


This process continues for the third, fourth and fifth level comparisons of all the 
ldrXt.out files and the same techniques are applied to all levels of comparison of the 


ldarXd.out files. 
2. Comparison of “T” and “D” Type Files 


The function m/d212 (for Mission 1 Duty 2 versus Treatment 2) was written to 
compare the most recent versions of the /dr2t.out (mlout2t) and ldr2d.out (m1out2d) 


files. Both files are read into the function and renamed /dr2t and /dr2d, respectively. A 
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sixth column is then added to each file. Since the patients on limited duty are unable to 
perform specific activities, the /dr2d.out is treated as the “main” file. As a result, all 
records in the /dr2d file are assigned a “1” in column six to indicate that particular patient 
is counted as “out” for one day. Each record for each day is then compared between the 
two files. If there is no match of records for a day, a “1” is placed in column six of the 
ldr2t file (so that it is counted for that day in the treatment file). If there is a match of 
records for a given day, a “O” is placed in column six, indicating it has been accounted for 
by the /dr2d file. Function m/d212 is given in Table 3.14. A sample output of the /dr21 file 
(named mJtrt2.out) after completing the steps described above is given in Table 3.15. 


function(m1 out2t, m1dout2d) 
{ 
# mission1: compare trmt Id2 to duty Id2 

# to check for equal entries 

# If the same person with the same diagnosis for 

# a given day is put on limited duty and is seeking 
# treatment, man-hrs lost will only be accounted for 


# in the limited duty file; if he's only being treated 


# man-hrs lost will be accounted for in the treatment 

# file 
vec.2d <- apply(m1dout2d[, c(1, 3, 4)], 1, paste, collapse = "”) 
vec.2t <- apply(m1out2t[, c(1, 3, 4)], 1, paste, collapse = "") 
matchers <- match(vec. 2t, vec.2d) 
Idr2d <- cbind(m1dout2d, rep("1", nrow(m1 dout2d))) 
Idr2t <- cbind(m1 out2t, rep("T’, nrow(m1 out2t))) 
# if it does match, put a 0 where the T is in the trmt file 
new.col <- ncol(idr2t) 
Idr2t{tis.na(matchers), new.col] <- "0" 
# if it does not match, put a 1 where the T is in the trmt file 
Idr2t[is.na(matchers), new.col] <- "1" 
return(list(idr2d = Idr2d, Idr2t = Idr2t)) 





Table 3.14 Function m/d2t2 to compare “treatment” and “duty” LD code 2 


Day _|Cat|Demo [Diagcode __|Ldcode [Cu 
[960212 |T__|204 _|724502000 2 


960220 |T_ (231 [72190000020 
960221 |T__|231__|721900000 1210 


Table 3.15 Sample output of mJ/¢rt2. out file 
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This comparison is made for the remaining LD codes 5, 6, 10, 11, and 12. The 
updated files are then assigned appropriate column headings (Day, Cat, Demo, Diagcode, 
Ctr). 


3. Computing the MOE 


The last step in computing the MOE is to perform the appropriate computations 
for each day using the updated files. Recall the assumptions for computing the man-days 


lost for patients receiving treatment and/or on limited duty: 


e Ifon limited duty only , multiply the “Ctr” entry by 1 if'a “1” appears ( entry 
of “1” for man “out” that day; entry of “0” for those accounted for by 
another LD code). 

e Ifreceiving treatment on/y , multiply the “Ctr” entry by 0.5 if a “1” appears 
(entry of “1” implies not accounted for by another LD code or by the 
comparable “duty” file; entry of “O” otherwise). 


The computations are performed by reading the updated files output by the 
aforementioned functions into another function to compute the day totals (i.e., function 
mIdaytotals computes the day totals for Mission 1). This function sums, across all tables 
relating to the mission or mission segment (e.g., tables generated for LD codes 2, 5, 6, 10, 
11, and 12 for Mission 1), the number of man-days lost each day due to individuals placed 
on limited duty and receiving treatment. Additional columns are created computing the 
man-hours lost per day due to placement on limited duty and receiving treatment. The last 
two columns created represent the total man-days lost and man-hours lost per day. The 
resulting file for Mission 1 (m/J/totals) is exported to Microsoft Excel and named 
mItotals.xls. All S-Plus commands used for computing the MOE for Mission | are given 
in Appendix D. A sample output of the file is given in Table 3.16. The complete file is 
given in Appendix E. 


Date__|D.Man-Days |D.Man-Hrs |T.Man-Days |T.Man-Hrs |Man-Days |Man-Hrs | 
Se0715|3 alla, |) el ec STS cco 


ge0ii6|4 tp Ss sd Sd 


Table 3.16 Sample output of m/totals.x/s file 
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4. Subsequent Computations 


As will be seen in Chapter IV, the graphical results of man-days lost over time for 
each mission and mission segment show little variation in totals from mission(s) to mission 
segment(s). This is primarily due to the fact that of the 73 unique diagnosis codes, 
approximately 44 percent could potentially result in LD code 5 (no upper body) and 
approximately 38 percent could potentially result in LD code 6 (no lower body). 
Although only seven percent of those codes are considered to result in LD code 2 (no 
duty), those diagnosis codes making up the seven percent occur frequently in the database 
provided. Asa result, these three LD codes have a significant impact on man-days lost. 

In order to visualize differences based on specific activities, the MOE was 
modified somewhat and new computations were made. The modified MOE is defined as 
activity dependent man-days lost over time due to MS injuries incurred by Navy SEALs. 
Three specific activities in which SEALs were limited in performance were chosen: 
diving, swimming with fins, and running/walking under stressful conditions (corresponding 
to LD codes 4, 9, and 12, respectively). 

This results in modification of the function which computes the day totals. The 
tables for LD categories 4, 9, and 12 are already existent in S-Plus. The corresponding 
“treatment” and “duty” tables were referred to as /drXt.out and ldrXd.out, where X 
represents 4, 9, and 12. Therefore, in computing man-days lost over time due to Activity 
4, only the /dr4t.out and /dr4d.out files have to be compared. The computations and 
creation of a table are performed by function daytotalsact4 (see Table 3.17). An additional 
column, column five, is added to each file. A “1” 1s placed in every entry of column five in 
the /dr4d. out file (considered to be the “main” file). If the same record (fields one through 
four) appears in both files, a “0” is placed in column five of the /dr4t.out file to indicate it 
was accounted for by the /dr4d.out file; otherwise, a “1” is placed there to indicate it 
should be accounted for by the /dr4t.out file. Sample output of the file generated by 
function daytotalsact4, a4totals (saved as a4totals.xls after export to Microsoft Excel), 
which is used to plot man-days lost over time due to Activity 4, 1s given in Table 3.18. 


The output file generated for Activity 4 is given in Appendix F. 
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a4output <- matrix(O, nrow = 436, ncol = 7) 
a4output[, 1] # now the first col is each of the 436 days in the file 
for(i in 1:nrow(a4output)) { 
day <-adout # first day is 951204 
tableda4 <- a4dupd[a4dupd[, 1] == day, ] 
dman.days <- sum(as.numeric(tableda4[, 5])) 


dman.hrs <- 24 * dman.days 
tableta4 <- a4tupd[a4tupd[, 1] == day, ] 


tman.days <- sum(as.numeric(tableta4[, 5])) 

tman.hrs <- 24 * tman.days 

a4output[i, 2] <- dman.days 

a4output[i, 3] <- dman.hrs 

a4output[i, 4] <- tman.days 

a4output[i, 5] <- tman.hrs 

a4output[i, 6] <- as numeric(a4output[i, 2]) + 0.5 * as.numeric(a4output[i, 4]) 

a4output[I, 7] <- as.numeric(a4output[i, 3]) + 0.5 * as. numeric(a4output[i, 5]) 
} 
return(list(a4output = a4output)) 





Table 3.17 Function daytotalsact4 to compute MOE as a function of “no diving”. 


[Day __|D.Man-Days |D.Man-Hrs |T.Man-Days |T.Man-Hrs |Man-Days |Man-Hrs | 
geozi2io fo SQ 4g 


960212] OI) | 0 6 aaa ae eee 
geozi4jo tots 8 


Table 3.18 Sample output of a4totals.x/s 





All graphical and statistical results for Mission 1, Mission 25, Mission segment 1, 
Mission segment 3, and Activities 4, 9, and 12, are presented in Chapter IV. 
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IV. RESULTS 


A. SAMPLE MISSIONS AND MISSION SEGMENTS 


High frequency graphs were produced for Missions 1 and 25 and Mission 
Segments | and 3 showing man-days lost over time, dependent on each of the missions 
and segments. As mentioned in Chapter III, the results show little variation across this 
sample, primarily due to the impact of LD codes 2, 5, and 6. Before displaying results, it is 
important to note which LD codes affect each of the missions and mission segments in the 


sample. This information is given in Table 4.1. 


LD Codes Affecting Each Type 


Table 4.1 LD Codes affecting each Mission/Segment in Sample 







Graphical results are presented in Figures 4.1 through 4.4. Statistical results for 
the mean, variance and standard deviation of the number of man-days lost per day due to 


each of the missions and segments are given in Table 4.2. 


| type | cet per bay” | Varianes [Deviation 
Type Lost Per Day Variance | Deviation 
| Mission25— | 12.11 3.35 | 56 


Table 4.2 Statistical Results - Missions and Mission Segments 







As can be seen in Table 4.1, there is little difference in the mean number of men 
“out” per day for each of the missions and mission segments in the sample. Of more 
significance is the fact that, on average, approximately 12 men were out each day 
regardless of which mission or mission segment was considered. Keeping in mind that 
there are approximately only 900 SEALs attached to NSWG-TWO, 343 of which visited 
one of the clinics for MS injuries throughout the 13 month period, this number may be 
significant. The facts that many of the 900 SEALs attached at NSWG-TWO are on 
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deployment six months over the course of a year and that probably some of the data were 


lost may account for only 343 actually receiving treatment at NSWG-TWO Medical. 


Mission 1 Man-days Lost Over Time 
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Figure 4.1 Mission 1 Man-days Lost Over Time 
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Mission 25 Man-days Lost Over Time 
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Figure 4.2 Mission 25 Man-days Lost Over Time 
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Mission Segment 1 Man-days Lost Over Time 
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Figure 4.3 Mission Segment 1 Man-days Lost Over Time 
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Mission Segment 3 Man-days Lost Over Time 
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Figure 4.4 Mission Segment 3 Man-days Lost Over Time 


B. SAMPLE ACTIVITIES 


In order to more clearly see the variation from day to day as a result of particular 
activities utilized either during training or as an integral part of performing a specific 
mission Or mission segment, the MOE was modified to be activity dependent man-days 
lost over time. Three different activities (diving, swimming with fins, and running/walking 
under stressful conditions), which are directly affected by LD codes 4, 9, and 12, 
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respectively, were chosen. Graphical results of each are displayed in Figures 4.5 through 


4.7. The statistical results are given in Table 4.3. 


acti Mean # Man-days 
Activit Lost Per Da Variance | Deviation 


Table 4.3 Srncical Results - Kenn 











It is evident, both graphically and statistically, that certain activities which SEALs 
perform regularly can lead to MS injuries and significantly impact the man-days lost over 
time as a result. 

For example, consider Mission | which 1s directly affected by LD codes 2, 5, 6, 10, 
11, and 12. As shown in Table 4.2, the mean number of man-days lost over time due to 
Mission | is 12.64. As indicated in Table 4.3, activity 12 (which 1s directly affected by LD 
code 12) resulted in 5.47 man-days lost over time. Of course, there would be some 
overlap by the other LD codes, but activity 12 alone makes up almost 44 percent of the 
average man-days lost over time for Mission 1. 

A similar explanation also holds when comparing the results for Activity 9 to 
Mission Segment | results. The mean number of man-days lost due to Activity 9 alone, 
given in Table 4.3, is 4.62, while the mean number of man-days lost due to Mission 
Segment 1, given in Table 4.2 is 11.49. This implies that 40 percent of all man-days lost 
due to Mission Segment 1 can be attributed to not being able to perform (or perform to 
maximum potential) Activity 9, swimming with fins. 

Finally, consider the comparison of results for Activity 4 and Mission Segment 3. 
As given in Table 4.3, 3.08 man-days lost are attributed to limited or no performance of 
Activity 4 and 12.38 man-days are lost due to Mission Segment 3 (Table 4.2). Because 
performance of Mission Segment 3 relies on being able to perform Activity 4, 
approximately 25 percent of the inability to perform Mission Segment 3 is a direct result 
of the inability to perform Activity 4, which 1s diving. 

Again, it is important to note that these are estimates, because the last three 
comparisons do not take into account the other LD codes associated with each of the 


missions and mission segments in the sample. 
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Activity 4 Man-days Lost Over Time 
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Figure 4.5 Activity 4 Man-days Lost Over Time 
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Activity 9 Man-days Lost Over Time 
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Figure 4.6 Activity 9 Man-days Lost Over Time 
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Activity 12 Man-days Lost Over Time 
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Figure 4.7 Activity 12 Man-days Lost Over Time 


C. OVERVIEW OF GRAPHICAL RESULTS 


1. Similarities in Overall Shape 


As seen in Figures 4.1 through 4.7, a majority of the man-days lost over time occur 
between days 90 and 190 (corresponds roughly to the period of the middle of March 
through the end of June). The majority of man-days lost during this time period may be 


43 


because there were more SEALs present at the command (i.e., not on deployment) during 
that time; and as a result, appeared more frequently in the database during this time. There 
are also peaks which occur at various times throughout the course of the 13 month period, 
which would also indicate some type of seasonal variation. 

There 1s little activity shown by any of the graphs at either the beginning or end of 
the time period of data collection. A possible reason for the small number of man-days 
lost per day in the beginning is the fact that the data collection started during the holiday 
season. Another reason may be that some data were lost and unrecovered. It was 
determined that the system of data collection, MADSYS, did “crash” on three separate 
occasions and that this occurred when the system was first put into use. 

The man-days lost over time tapers off toward the end of the period of data 
collection and coincides with the fact there was little activity at the same time in 1996. It 
is unlikely that any data were lost during the last few months of data collection, because 
the person at NSWG-TWO who input the data was careful to save the information more 
regularly and also because the system had not “crashed” since April, 1996. 


2. Differences in Overall Shape 


Both Activities 4 (diving) and 9 (swimming with fins) are aquatic activities. 
However, there appears to be a peak in the graph of Activity 4 man-days lost over time 
(Figure 4.5) from day 260 through 305 (corresponding to mid June through August of 
1996), which does not appear in the graph of Activity 9 man-days lost over time. This 
may be due to the fact that diving relies more on the use of the lower body, primarily the 
use of the hips and thighs. However, surface swimming involves greater reliance on the 
upper body, thus reducing the demands on the lower body. 

Types of data collection necessary to more accurately determine man-days lost 
over time (due to missions, mission segments, or specific activities) are discussed in the 


next chapter. Also, areas of recommended future research are described. 
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V. RECOMMENDATIONS AND CONCLUSIONS 


A. SUGGESTED DATABASE DESIGN 


1. Data Items to be Collected 


In addition to the data items collected by the software system MADSYS, other 
items should be collected. As mentioned in Chapter III, several data items were missing 
from the database used for analysis in this thesis. For example, the healing time estimate 
for a particular patient with a specific injury was not indicated. This information would be 
particularly useful if one had a computer program which would allow the user to specify a 
day, and output not only the current man-days lost, but also predict the number of man- 
days lost for some time in the future - based on the healing time estimate recorded by the 
reviewing physician. 

Furthermore, there should be an indication of the patient’s date of full recovery for 
a particular diagnosis. Ifa patient is still injured at the time of transfer to another location, 
separation from the service, or any other reason, this information should be recorded as 
well. It may also be useful to indicate when an individual is preparing for deployment in 
order to continue tracking and updating patient information. 

In order to obtain a realistic report of man-days lost over time due to injury, the 
reviewing physician, together with the patient, should include which types of activities the 
individual either would not be able to perform at all or would result in a reduced 
performance capability. This information is especially important because certain injuries 
affect individuals in different ways as mentioned in Chapter III. One person’s injury may 
result in his ability to perform a particular activity more effectively than another individual 
with the same injury. Activities considered essential that SEALs be able to perform 
should be listed and given a corresponding percentage of performance ability based on an 
agreement between the reviewing physician and the patient. Those activities given a 
“zero” percentage of performance ability would, obviously, refer to those in which the 
patient received a limited duty status by the physician. Those given a percentage other 
than zero would indicate that the individual would still be able to perform the activity, but 
not to his maximum capability. As time progresses, with the patient receiving treatment, 


this percentage would likely increase for a particular injury. Therefore, in order to 
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compute man-days lost more accurately for any given day, this information should be 
updated - at least with every visit in which a patient receives treatment for that diagnosis. 

Therefore, rather than creating a mapping from diagnosis code to LD code 
regardless of the individual involved (as was the case in this thesis), this method of data 
collection would show how each individual 1s affected by each injury. This would allow 
for further analysis to determine if patients with the same diagnosis were affected in the 
same manner. 

Another important consideration is the fact that many patients in the database 
were not given limited duty by the reviewing physician, but were receiving treatment 
frequently for long periods of time. Therefore, performing any activities which would 
directly utilize the muscles and joints injured would put the individual at risk. Determining 
the risk factors involved is another large project itself, but should be given consideration. 
Obviously, if an injured person continues to perform activities which directly affect the 
muscles and joints injured, this could prolong the injury and may even result in limited 
ability to heal properly. This could then lead to recurring problems and further complicate 
the matter. Because the SEALs are required to perform the most physically demanding 
types of activities, this population is even more at risk. If the Navy wants to keep those 
individuals who are so hard to acquire in the first place, and to keep them performing to 
their maximum potential, every precaution should be taken to minimize the risk of injury 
and to provide the best care in order that the affected areas heal properly and in a timely 
manner. 

This suggestion further supports other work Dr. Prusaczyk had been conducting in 
collaboration with NSW medical personnel: to install facilities such as sports clinics, 
physical therapy clinics, and clinics housed with athletic trainers in NSW commands where 
such facilities do not exist. Not only will it be beneficial to the Navy in collecting the data 
items necessary to properly track man-days lost due to MS injuries incurred, but the Navy 
could also benefit by providing better care for its sailors whenever possible. 

This type of data collection is not limited to Navy SEALs with MS injuries. The 
methodology can be implemented for all military personnel sustaining any type of injury 
which, as mentioned earlier, is currently lacking. To further stress the importance to the 
Department of Defense (DoD), adding the rank and associated pay scale for each patient 
in the database, along with the other parameters used in computing the MOE, would be 
beneficial. This would allow the users of the system to see on any particular day the actual 


cost to the DoD of personnel out with injuries that day. 
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2. Use of a Central Database 


A centralized system would retrieve information on a daily basis from all Medical 
Facilities under its cognizance. The central location, of course, would be the controlling 
center. Each command under it would use the same system of data collection which 
would “be fed” into the central computer. In addition to the controlling center being able 
to determine aggregated operational information over all of its affiliate commands, the 


affiliate commands would also be able to assess their operational status at any time. 


B. RECOMMENDATIONS FOR FUTURE RESEARCH 


1. Software Design 


Obviously, the design of the database as mentioned above constitutes a very 
important and practical topic for future research. In addition to actually writing a program 
to develop the database, there is a need to automatically update the MOE, perhaps based 
on a date input by the user. Furthermore, a program which would also forecast, based on 
the healing time estimates, the man-days lost for cost to DoD from a given date is needed. 
The entire program could encompass both aspects: generating the data on spreadsheets 
easily accessible by the user and generating a “cost to DoD” estimate either for the current 
time or some future time period. 

With a system as described above, medical professionals would be equipped to 
determine the usefulness and effectiveness of their facilities in terms of providing health 
care for patients. Such information would aid in the ability of the facilities to determine 
areas where improvement is needed. The entire career of a SEAL could be tracked and 
researchers could determine if there are any injury patterns which individuals followed 
throughout the years. For example, if a patient encountered an Anterior Cruciate 
Ligament (ACL) tear (a particular type of knee injury) early in his career, the likelihood of 
experiencing similar injuries in the future could be established. Answering questions of 
this nature would greatly aid medical professionals in reassessing how to improve their 


role in the individual health care process. 
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2. Determining Distribution of MS Injuries 


Once the proper data items were collected using the new database design, fitting 
the man-days lost over time to a particular distribution would allow researchers to 
determine if any seasonal variation exists and provide them with the proper statistical 
means for forecasting based on the particular distribution. Many other types of analyses 
could be conducted such as determining if there exists a relationship between the activities 
which resulted in particular injuries and the injuries themselves, or if individuals sustaining 
the same type of injury were affected in the same manner. If the data revealed both the 
physician’s and individual’s assessment of percentage of performance ability, analysis 
could be conducted to determine if there exists any correlation between the two. 

The frequency of occurrence of injuries could provide guidance as to which 
injuries need to be more closely examined in terms of providing care. If researchers know 
which injuries are most likely to occur, specific action can be taken to reduce the 
occurrence of such injuries. For example, if the most frequent injuries are back-related 
injuries, then additional precautions may be taken at the training level to insure that 
proper exercise and strengthening of the back occur more frequently to minimize the 
occurrence of the injury. For those individuals in an operational status, the same 
techniques may be implemented to reduce occurrence of injury. 

Finally, a ratio MOE could be computed for each affiliate command by dividing the 
man-days lost for each day by the total man-days available for any specified mission, 
mission segment or activity. This measure would provide a dynamic indicator of the 


operational readiness of the command to perform its assigned missions. 
C. CONCLUSIONS 


As stated earlier, MS injuries have been one of the major areas of research at 
NHRC, within the NSW community, and also within DoD. The creation of a central 
database and affiliate command databases for the tracking of MS injuries will provide the 
DoD with several important pieces of information. Such information includes how to 
improve each of its medical facilities based on an assessment of the effectiveness of each 
facility. The information retrieved from the database will also provide a means of 
incorporating better training and exercises to minimize risk to personnel. The information 


provided will also be useful as a contingency planning tool by recording the impact on 
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operational readiness, both in terms of man-days lost over time and estimation of cost to 


DoD, as a result of injury incurred by military personnel. 
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APPENDIX A. SAMPLE OUTPUT OF TREATMENT FILE CREATED IN 














PARADOX 
INOXD|NOXV|NOXD|NOXT|TRIMJENCDATE [CLINI |VISTY|ICOT [INJTY[INyme[AcTc [acts [Sex [RACE [BRANCH | 
|_102] 3829] 981| S402] sjivsoe_ [ttt 
|_102] 3830] 981] 5405] sfiv7ee_ | at tt tt 
| 102] 3830] 981] S404] 6tivee_ | tt 
|_102] 3830} 981] S403} tives | ttt 
| 105] 105f 105] tor zee | ttt 
| 105] 110] 105] 105] ofsye_ | ttt tt tt 
| 105] 110} 105] 103] iss | ttt tt 
| 105] 110] 105] 104 ates | ttt tt 
| 105] 299] 264] 24a} ative | att ttt 
| 105| 299] 264] 243|s{vivee | att tt 
325] oizes | at tt 
| 105] 557] 264) 325] afavigvo6_ =| Tt tt 
| 105] 574] 105] 326] gja/iz96_ | ttt 


| 105] 574{ 105} 326] 92/1396 | St ttt 
| 105] 575] 105] 327] gjaviaop_— | at ttt 
| 105] 575{ 105] _327|_afziaoe_— | at tt tt 
| 105] 576] 105] 328] ajanee_ | a] ttt 
| 105] 576] 105] 328] afareoe | at tt ttt tat 
| 105] ses] 105] 329] 9396 =| at ttf Stat 











DATEDEMO [{ENCDATE.1 |DIAGCODE |DIAGDESC DIAGSPEC DIAG DATE 










17296 —*|11/596_| 726691000 
1/2/96 111/596 __|726691000 
1/2196 2/22/96 
1/2/96 1/2/96 
17296 «| 1/3/96____—| 844200000 
8 
1/2/96 


Note: The two tables combined by having the lower table directly to the right of the upper 
table represents the first 20 records of the treatment.xls file. Four columns (representing 
last name, first name, middle initial, and social security number) have been eliminated for 
privacy considerations. Some column headings are hidden in order to be able to fit the 
table on the page. The first encounter date (ENCDATE) is the date the patient visited a 
particular clinic for an injury. The second encounter date (ENCDATE. 1) is the date of the 
first visit a patient had for a particular diagnosis and was not used in any of the analysis. 
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APPENDIX B. TRTDAYS AND LIMDUTY FILES CREATED IN S-PLUS 
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231 ssdt4a3 sdf 25g 281 dG SCC‘d‘ 15960315 {960329 1721900000 
273, 4491315 285 CC‘L'100._:—«#( 960306 __ 960315 |845001B00__| 


135 dt 554 36 Sid 2875 Cid 931960319 1960619 _|84080A000___ 
144598 147d gg 2 C*d 21 (960319 1960408 _|72889D000__ 
292 sdf 628d 335 Sid 292 Ss? SC“ Cid IS [960318 9604018) 17266970 R 
294631 337 Sid 294 Cd 7960319 1960325 1724502000 ___ 
298641344296 dd 8960320 1960327 __|728710000___ 
313,678 sdf 363_ 305 5 Cd 8960325 1960401 1840000000 
314s 680s 364.306 S77 Cid 8960325 1960401 __ 1728710000 ___ 


327. ~~ ~=|702 +=+{378 1316 #$.|j8 | 6 {|960327 726120000 
328 1705 |380 %|{317  #+x,|[5 {17 |960313 [960329 [727050000 
329 ~~ we Sov [3816 el 208e es |S ae  I.81960322 


329 [712 +=|381 +=+%|319 {3 {9 {960328 |960405 [840000000 _ | 
343. {906 |[398 [324 #+x,|5  —_|9 [960328 [960405 |922000000_ | 
344 —-f[951 ~—s [399 ~S 325 Ss f0—S—<$s sd 16 ~=1960403 (960418 |84200B000_ | 
345 1953 [400 [325 [5 [16 [960403 [960418 {72340A000 _ 
294. ~«4(1019~—=s«4} 337, =Sssdf340)—0—dT 3 Sid 13“ (960417 1960429 [724502000 | 
357 1025 420 1342, 16~—~é<“—sS Ss 4960418 1960501 _ 1 726691000 

184. —=—«({1050_ sis (345. 16 —*«| 6 [960411 _|960416 |726710000 _ 
362, 11108 [427 [348 [6 [14 [960403 |960416 {844200000 _ _ 
368.1155 1434 (351 S38 sd 15 1960417 1960501 
369 (1158 [435 [352 [7 [152 |960208 |960708 {844200000 

373. 11208 [440 [355 |7 [31 [960423 {960523 |845001B00_| 
378. 11245 [445 [360 I[5 [14 [960423 [960506 |726191000__| 
385 11282 {454 |364 [5 __|32 {960311 _ |960411 
385 11283 [454 |366 |{5 [18 {960426 
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703 1960801 _|172889A000 
960724 |845001A00 
1836000000 


557 16st 8 (1960814 1960821 {836100000 __| 
ip59_—s*d'2611 Ss 733571 Sid C*dT' 1,:—«SYS' 960814 _[960824 {922000000 ___ 
572_—Ss—=«d(2725”—s 760 Sst 598 =| Cid 11.:=—«4' 960827 1960906 __|84080A000___ 
574 s—«d(2738~—sdt 762s 600, =F C“Ci;i‘L I= (960826 {960910 |845001B00__| 
s87_——s*d'2834_—st 7931632 Cd 93 1960820 961120 __|726191000__| 
641__s—«d}3369 Ss 898_ sf 739 Sid 7 SC —Cid 11 —« (961008 __ 1961018 _|82380A000__ 
643, ——s*d('3396 902, sf 74a_ tS Cd 23, [1961010 __ 1961101 __|84212A000___ 
647_s—sdi3a59_—s 909s 750 Cid 7_—«('961016 __ 1961022 1727100000 __| 
6603569. sf 930s 768 CE 9610231961031 _171946A000__| 
666 —«s(3631_~—sf 943,780 7 id 16 [961030 __]961114 1844201000 __| 
683 s«d(3729| Sst 980s 816 5 15 1961108 1961122 1831000000 __| 
697_—«d3939_—sdf1006_—_—“_ at s43 Sd 26 [961118 1961213 1811000000 __| 
6853983987 st 851 dd A f'967122 [961125 1923200000 ___ 
705 s«s{4o70_—sdf1026~Ss—idt 863 CC‘ 7-~=—=«d('961125__ 1961201 1836000000 __| 
6714081951 sf s70_—S id Cid 1881961126 __ 1961213 _|724502000___ 
716___—s«4j4201__—sdt1050_~sdts89_ 7] Cid 151961205 __|961219 1845001800 _| 
731 —idja240_—sdtn080._~ssdt 922, 7] Cid 981961210 _]970317__|844200000 
765 s*i4701_~—sidw176~Ssidtn018 = s*d'196 [970107 _]970721 1924000000 
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APPENDIX C. FILE REFLECTING MAPPING OF DIAGNOSIS CODES TO LD 
CODES 


Una a Vs 
al ble alee 1 a 
715970000] _ o| 1| o| o| o ei 
719410000] “o|ro| w1| “ol ol o| {a ati ol ol 
/71946A000} o| i| of of of of of of of of 14 
| 721900000} 1| of of of 1/ of of of of of of 
| 722000000 _o| 1] of 11 1/ of of of of oo] of 
722100000] o| 1| of of 1] of of of 1 of oO 
|722400000] 0} o| o| 1| o}| of o| of a1] of of 
| 723100000} _o| of o| 14| of} o| o| of 4] o| of 
| 723400000] o| ol o| 1| o}| of of} of 4] o| of 
|72340A000] o| o| o| 1| of} o| 1| 4] 14/ o| oo 
| 724501000] o| o| o| o| 14] o| of} o| of} o| oj 
| 724502000] 1 1| of 1] 1] of of of oo] o| o 
| 726000000] _o} o| 14] 1| of of 14] 14] 14) 4] Oo 
|726100000} _o| o| 1] 1/ o| o| of of 4) 4] OF 
726120000] o| o}| o} 1| of of 14] a1/ 14] of oo 
| 726191000} o| o| 14] 1) o| of o| of 4] of of 
| 726310000/_o| o} o| of o| of} of o| a4| o| o 
| 726320000] o| o} o| 1] o| o| of} o| of of o 
| 726390000|_ 0| oO} o| o| o| of} of o| 1| o| of 
| 72664A000{_ o| o| o| of} 1] 14] o| o| oo] of} 11 
| 726650000/_ o| o} of o| a] 14] 1) 4] of o| 11 
| 726691000] o| o| o| of 14] 1] of 1] of of} 14 
| 726692000] _o| of; o} 1] of of 1] 14] 4] of oO 
| 726740000 ewio eto Mei} ewig | ai rel 1G | NG a 
726799000 6) Go} oll sol *o| a ol “Po a 
| 727050000] o/ o| o| 1 of} o| o| of of of of 
| 727100000] o| o| of o| o| o| of} o| o| o| of 
| 728710000] o} o}| of} of 11 41 o| of oo] o| 11 
72889000] _o| o| 1] 4/ of} o| 4] «| a4] 4| OO 
-72889D000/_1| o| o| of of o| of o| of o| o 
739100000] o| o| o| o| of of of of of of o 
739200000] o| o| o| o| oj} of o| of of of o 
739300000] _o| o} o| o| o} o} o| of o| of o 
-805soo000] _1| o| o| 1] of of o| o| of o| o 
811000000] _o| o| of 1| o| o| a| 1] 4] of oo 
813800000] o| o| of «| o| of o| of 4] o| oO 
813900000! ol o}| of 1] o}| of} «| 1) 4) o| oO 
|g1so00000| of o| o| 11 o} of of of 4] of} 
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gioco) ol ol ilo Ge ae ood oleae 
| 822000000] 0} o| of of 1] a] a] af of of 74 
| 823800000] of of of of 1} 1] of of of of 4 
| 82380A000] 0] o| o| o| o| 4] of of of of 141 
823900000] of o| of o| 1| a1 of of of of 14 
| 824800000] o| of o| of 1] 4] 14] 4] of of 1! 
| 825290000 0} o| of o} 1] of of 1] of of 4 
| 831000000] of 1| 4] 14] of of of of 4] Oo 

| 83100000] of of 1| 4] o| of o| of 41 of oO 
| 832000000] of o| of 1{| of of 1] 4] 1] of oO 
| 836000000] o| 14{ of of 14] of of 4] of of 14 
| 836100000] of 1| of of 14] 4] of of of of 14 
|840000000]_o| of 1| 4] of of 4 4] 1) of Oo 
/84080aooo| _o| of 1| 4/ o} of} of of 4} of oO 
| 840808000] _o} of 1| of o| of of of 4] o| oOo 
84190000] o| o| o| 4] o| oj} o} of 4] o| 
| 842008000] _o| o| o| 1] of of of of} 4) of oo 
| 84212A000| 0] o| o| 14 o| of} of of 11 of o 
| 843808000| 0} o| o| o| 1] of of of of 4] Oo 
| 844000000] _o} o| o| of 1] of of of of of 14 
| 844100000] 0} 1| of oo} 1] 4] of of of of 4 
| 8442000000} o| of of 1] 4] of 1] of of 14 
| 844201000|_ 0} o| of o| of 1 o| o| of of 14 
|84480A000| ol of of of 1| 4] o| of oo} oo} 14 
| 845001A00/ 0} o| of o| of 1] of 1] of} of 14 
/845001B00| of 1] of of 4] a{ o| of o| o| 14 
847100000] _o} o} of} 1{ of of of of 4] 4] Oo 
|84aaog000| of of of of 1] a{ 4] 4] of oo} 14 
|922000000/ _o| of of 1] of o| of of 4] 4] Oo 
|-923100000] of o| of a] of of o| of 4] of oO 

a) 






























62 


APPENDIX D. S-PLUS COMMANDS AND FUNCTIONS FOR COMPUTING 
MISSION ONE MOE 


Mission 1: 


># Requires activities affected by LD codes 2, 5, 6, 10, 11, 12 

># Will makes pairwise comparisons for all treatment file information 

># which are (2,5), (2,6), (2,10), (2,11), (2,12), (5,6), (5,10), (5,11), (5,12), 

># (6,10), (6,11), (6,12), (10,11), (10,12), (11,12) 

># This will require 15 steps for comparison of treatment file information 

># Function names are m1tl, m1t2, .... m1t15 referring to Mission 1 Treatment Step 1,... Mission 1 
># Treatment Step 15 


># Creating the Tables IdrX.out, where X represents 2,5,6,10,11,12 
> # First look at the function which will create a table for LD Code 2; 
> # Similar functions are used for creating tables for the other LD Codes 
> Idtable2 
function(rowcodes, unqcode) 
{ 
# rowcodes here is the dutyntrt file 
# unqcode is the unqcode file 
t2.ctr <- 1 
table2 <- matnix(" '", nrow = 2500, ncol = 4) 
rowcodes.col <- ncol(rowcodes) 
for(i in 1:nrow(rowcodes)) { 
cat("Now on row ", 1) 
j.seq <- seq(4, rowcodes.col, by = 3) 
jS <- rowcodes[i, j.seq] 
js <- jsfjs != "99999" & js !=" 99999"| 
j.seq <- seq(4, by = 3, length = length(s)) 
samprows <- unqcode[as.logical(match(unqcode[, 1], js, 0)), , drop = F] 
which.are.2 <- samprows[samprows|, 2] == "1", 1] 
if(length(which.are.2) != 0) { 
cat("...found!\n") 
j.seq <- j.seq[as.logical(match(js, which.are.2, 0))] 
right <- cbind(i, c(rbind(1, j.seq - 2, j.seq - 1, j.seq))) 
n <- length(j.seq) 
left <- cbind(rep(t2.ctr:(t2.ctr + n - 1), rep(4, n)), rep(1:4, n)) 
table2[left] <- rowcodes[right] 
t2.ctr <- t2.ctr+n 


} 
else cat("...nothing there\n") 


} 
return(list(table2 = table2)) 


j 
> Idr2out_Idtable2(dutyntrt,unqcode) # passes in the files dutyntrt and unqcode 
>Idr2.out_Idr2out$table2 # extracts the table 2 created by the function, stores it as Idr2.out 


>Idr2t.out_ldr2.out[Idr2.out[,2|=="T”,] # extracts only the treatment records from the table 
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># The same commands are applied to each of the functions which are used to create tables 
># for the other LD Codes. The result is 5 files named IdrSt.out, Idr6t.out, Idr10t. out, 
># Idrilt.out, ldr12t.out 


># Now, the first level of pairwise comparisons , (2,5), (2,6), (2,10), (2,11), (2,12): 
># First, to compare (2,5) 
># Using the function ml1t1 


>mitl 
function(Idr2t.out, IdrSt.out) 
{ 
# 
# mission | tl: deluxe version 
# 
vec.2 <- apply(Idr2t.out, 1, paste, collapse = "") 
vec.5 <- apply(ldr5t.out, 1, paste, collapse = "") 
matchers <- match(vec.5, vec.2) 
Idr2t <- cbind(Idr2t.out, rep("2", nrow(idr2t.out))) 
Idr5t <- cbind(Idr5t.out, rep("5", nrow(IdrSt.out))) 
new.col <- ncol(Idr5t) 
Idr5t[{!is.na(matchers), new.col] <- "2" 
return(list(Idr2t = Idr2t, ldrSt = Idr5t)) 
} 
>m1.out_m1tl(Idr2t.out,ldr5t.out) 
>mlout2t_m1].out$ldr2t # Stores the resulting Idr2t table in mlout2t 
>mlout5t_m1.out$ldr5t # stores the resulting IdrSt table in mloutSt 
>m2.out_m1t2(Idr2t.out,ldr6t.out) 
>m2out6t_m2.out$ldrét # stores the resulting Idr6t table in m2out6t 


> # no need to keep outputting the ldr2t table because it will not change 
>m3.out_m1t3(ldr2t.out,ldr10t.out) 

>m3out10t_m3.out$ldr10t 

>m-4.out_m1t4(idr2t.out,ldr1 1t.out) 

>md4out] 1t_m4.out$ldr1 It 

>m5.out_m1t5(Idr2t.out,ldr12t.out) 

>mS5outl2t_m5.out$ldr12t 


># Second level pairwise comparisons: (5,6), (5,10), (5,11), (5,12) 

># Using the updated tables mloutSt, m2out6t, m3out10t, m4out11t, mSout12t 

> ml1t6 

function(m1out5t, m2out6t) 

{ 

4 

# mission | t6: deluxe version 

# This function receives the first updated version of 

# the file IdrS5t.out (mlout5t) and uses it to compare to 

# most recent version of Idr6t.out (m2out6t) 

# Can only use apply on cols | thru 4 of mlout5t and m2out6t 

# Places a "5" in col 5 of m2out6t if that row of m2out6t 

# exactly matches a row in mlout5t 
vec.5 <- apply(mlout5t[, c(1:4)], 1, paste, collapse = "") 
vec.6 <- apply(m2out6t[, c(1:4)], 1, paste, collapse = "") 
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matchers <- match(vec.6, vec.5) 
Idr6t <- m2out6t 
new.col <- ncol(ldré6t) 
not2s <- Idrét{!is.na(matchers),new.col] != 2 
# take the entries that haven’t already been accounted 
# for in ld2 file and assign them a “5” in col 5 
Idr6t{(!is.na(matchers)) & (not2s), new.col] <- "5" 
return(list(ldr6t = Idr6t)) 
j 
>m6.out_m1t6(mlout5t,m2out6t) 
>mé6out6t_m6.out$ldrét 
>m7.out_m1t7(mlout5t,m3out10t) 
>m7outl0t_m7.out$ldr10t 
>m8.out_m1t8(mlout5t,m4out1 It) 
>m8out1 1t_m8.out$ldr1 It 
>m9.out_m1t9(mloutSt,mSout12t) 


># Third level pairwise Comparisons: (6,10), (6,11), (6,12) 
># Using updated tables m6out6t, m7out10t, m8out11t, m9out12t 
>m1tlO 
function(mé6out6t, m7out10t) 
{ 
# 
# mission | tl0: deluxe version 
# This function receives the 2nd updated version of 
# the file ldr6t.out (m6out6t) and uses it to compare to 
# most recent version of Idr10t.out (m7out10t) 
# Can only use apply on cols | thru 4 of m6outét and m7out10t 
# Places a "6" in col 5 of m7out10t if that row of m7out10t 
# exactly matches a row in m6out6t 
vec.6 <- apply(m6outét[, c(1:4)], 1, paste, collapse = "") 
vec. 10 <- apply(m7out 1Ot[, c(1:4)], 1, paste, collapse = "") 
matchers <- match(vec. 10, vec.6) 
Idr10t <- m7out10t 
new.col <- ncol(Idr10t) 
not2s <- Ildr1l0t{!is.na(matchers),new.col] != 2 
not5s <- Idr10t{!is.na(matchers),new.col] != 5 
# take the entries that have already been accounted for by 
# 1d2 or Id5 file and assign them a “6” in col 5 
Idr10t{(!is.na(matchers)) & (not2s) & (not5s), new.col} <- "6" 
return(list(idr10t = Idr10t)) 
} 
>m10.out_m1tl0(m6out6t,m7out10t) 
>m10out10t_m10.out$ldr10t 
>m11.out_m1tl1(m6outét,m8out! It) 
>m1loutl1t_m11.out$ldr11t 
>m12.out_m1t12(m6out6t,m9out12t) 
>m12outl2t_m12.out$ldr12t 


# Fourth level pairwise comparisons: (10,11), (11,12) 
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# Using updated tables m10out10t, m11lout11t, m12out12t 
>m1tl3 
function(m10out10t, m1 lout] 1t) 
{ 
# 
# mission | t13: deluxe version 
# This function receives the 3rd updated version of 
# the file Idr10t.out (m10out10t) and uses it to compare to 
# most recent version of Idr11t.out (m1 lout] It) 
# Can only use apply on cols | thru 4 of m10out10t and m1 lout1 It 
# Places a "10" in col 5 of ml lout1 It if that row of m1 lout1 It 
# exactly matches a row in m10out10t 
vec. 10 <- apply(m10out 10t[, c(1:4)], 1, paste, collapse = "") 
vec.11 <- apply(m1 lout! 1t[, c(1:4)], 1, paste, collapse = "") 
matchers <- match(vec.11, vec.10) 
Idr11t <- mlloutllt 
new.col <- ncol(idr1 1t) 
not2s <- Idr1 It[!1s.na(matchers),new.col] != 2 
not5s <- Idrl 1t[!is.na(matchers),new.col] != 5 
not6s <- Idr1 1t[{!is.na(matchers),new.col] != 6 
# take entries that haven’t already been accounted for 
# by 1d2, Id5 or Id6 files and assign them a “10” in col 5 
Idr1 1t{(!1s.na(matchers)) & (not2s) & (not5s) & (not6s), new.col]} <- "10" 
return(list(idrl 1t = Idr11t)) 
} 
>m13.out_m1t13(m10out10t,m1 lout! 1t) 
>m13outl1t_m13.out$ldr1 lt 
>m14.out_m1t!4(m10out10t,m12o0ut12t) 
>m14outl2t_m14.out$ldr12t 


# Last comparison: (11,12) 
# Using updated tables m13out11t, m14out12t 
> mltl5 
function(m13out1 1t, m1l4out]2t) 
{ 
# 
# mission 1 t15: deluxe version 
# This function receives the 4th updated version of 
# the file Idrl1t.out (m13out11t) and uses it to compare to 
# most recent version of Idr12t.out (m14out12t) 
# Can only use apply on cols 1 thru 4 of m13out1 1t and m1]4out12t 
# Places a"11" in col 5 of m]l4out12t if that row of ml4out12t 
# exactly matches a row in ml3out1 It 
vec.11 <- apply(m13out1 1t[, c(1:4)], 1, paste, collapse = "") 
vec.12 <- apply(m14out12t{, c(1:4)], 1, paste, collapse = "") 
matchers <- match(vec.12, vec. 11) 
Idr12t <- ml4out12t 
new.col <- ncol(ldr12t) 
not2s <- Idr12t[!is.na(matchers),new.col] != 2 
not5s <- Idr12t[!is.na(matchers),new.col] != 5 
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not6s <- Idr12t[!is.na(matchers) new.col] != 6 
notlOs <- Idr]2t{!is.na(matchers),new.col] != 10 
# take entries that haven’t already been accounted for 
# by Id2, 1d5, 1d6 or 1d10 files and assign them an “11” in col 5 
Idr12t[!is.na(matchers), new.col] <- "11" 
return(list(ldr12t = Idr12t)) 
j 
>m15.out_m1tl5¢m13out] 1t,m14out12t) 
>m15out]2t_m15.out$ldr12t 


># Updated “treatment” tables are mlout2t, mloutSt, m6outé6t, m10out10t, m13out11t, m1S5out12¢t 


># Use similar functions for comparing all the IdrXd.out files 
># Results are the updated “duty” tables midout2d, m1doutSd, m6dout6d, m10dout10d, 
># mi3doutl1ld, m1Sdout12d 


># Now compare the updated “treatment” and “duty” tables 

># Function m1d2t2 compares mlout2t and midout2d: 

> m1ld2t2 

function(m]out2t, mldout2d) 

{ 

# mission1: compare trmt 1d2 to duty Id2 

# to check for equal entries 

# If the same person with the same diagnosis for 

# a given day is put on limited duty and is seeking 

# treatment, man-hrs lost will only be accounted for 

# in the limited duty file; if he's only being treated 

# man-hrs lost will be accounted for in the treatment 

# file 
vec.2d <- apply(ml1dout2d[, c(1, 3, 4)], 1, paste, collapse = "") 
vec.2t <- apply(mlout2t[, c(1, 3, 4)], 1, paste, collapse = "") 
matchers <- match(vec.2t, vec.2d) 
Idr2d <- cbind(mIdout2d, rep("1", nrow(m1dout2d))) 
Idr2t <- cbind(mlout2t, rep("T", nrow(mlout2t))) 
# if it does match, put a 0 where the T is in the trmt file 
new.col <- ncol(Idr2t) 
Idr2t[!is.na(matchers), new.col] <- "0" 
# if it does not match, put a 1 where the T is in the trmt file 
Idr2t[is.na(matchers), new.col] <- "1" 
return(list(ldr2d = Idr2d, Idr2t = Idr2t)) 

} 

>m1d2t2.out_m1d2t2(m1out2t,m1dout2d) 

>m1duty2.out_m1d2t2.out$ldr2d 

>ml1trt2.out_m1d2t2.out$ldr2t 


># Similar functions and commands used to obtain the following tables 

># midutyS.out, mltrtS.out, mlduty6.out, m1ltrt6.out, miduty10.out, mitrt10.out, 
># midutyll.out, mltrt11.out, mlduty12.out, and m1trt12.out 

># All these files are exported to Microsoft Excel. The filtering tool is used there 
># to determine which days in which SEALs were out should be counted by 
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># each one of the tables. 

># These updated files are then read back into S-Plus 

># Resulting files are named m1d2upd, m1t2upd, m1 dSupd, m1tSupd, m1d6upd, m1t6upd, 
># midl0upd, mitl0upd, midilupd, mitilupd, mid12upd, mitl2upd 


># Now reading to compute the day totals for Mission 1 using function daytotalsm1 
> daytotalsm1 
function?) 
{ 
mloutput <- matrix(0, nrow = 436, ncol = 7) 
mloutput[, 1] <- dutyntrt[, 1] # now the first col is each of the 436 days in the file 
for(i in 1:nrow(mloutput)) { 
day <- mloutput[1, 1] # first day is 951204 
tabled2 <- mld2upd[m1d2upd[, 1] == day, ] 
tabled5 <- mldSupd{[m1d5upd], 1] == day, | 
tabled6 <- mld6upd[m1d6upd[, 1] == day, ] 
tabled12 <- mld12upd[m1d12upd[, 1] == day, |] 
tabledm1 <- rbind(tabled2, tabledS, tabled6, tabled12) 
dman.days <- sum(as.numeric(tabledm1[, 7])) 
dman.hrs <- sum(as.numeric(tabledm 1[, 8])) 
tablet2 <- mlt2upd[m1t2upd[, 1] == day, ] 
tablet5 <- mltS5upd[m1tSupd[, 1] == day, |] 
tablet6 <- mlt6upd[m1t6upd[, 1] == day, ] 
tablet10 <- mltlOupd[{[m1t10Oupd[, 1] == day, ] 
tablet12 <- mltl2upd[m1t1l2upd[, 1] == day, | 
tabletm1 <- rbind(tablet2, tablet5, tablet6, tablet10, tablet12) 
tman.days <- sum(as.numeric(tabletm1[, 7])) 
tman.hrs <- sum(as.numeric(tabletm1[, 8])) 
mloutput{i, 2] <- dman.days 
mloutput{i, 3] <- dman.hrs 
mloutput[i, +] <- tman.days 
mloutput[i, 5] <- tman.hrs 
mloutput[i, 6] <- as. numeric(mloutput[i, 2]) + as.numeric(m1]output([1, 4]) 
mloutput[i, 7] <- as.numeric(mloutput[1, 3]) + as.numeric(mloutput[i, 5]) 
} 
return(list(mloutput = mloutput)) 
} 
>mIldaytotals_daytotalsm1() 
>mItotals_m1daytotals$m1output 


># Now plot Mission 1 man-days lost over time 

># First, create a vector of 436 days corresponding to dates 951204 through 970211 
>mIdates_seq(1,436,1) 

>win.graph() 

>plot(as.numeric(m|dates), as.numeric(m1totals[,6]), type = “h”, xlab= “Days 1 through 436, 
>continue string: (representing dates 951204 - 970211)”, ylab= “Man-days Lost’) 

> title(“Mission | Man-days Lost Over Time”) 
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APPENDIX E. MISSION ONE FILE: MAN-DAYS LOST OVER TIME 
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APPENDIX F. ACTIVITY FOUR FILE: MAN-DAYS LOST OVER TIME 
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